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Abstract 


Muon g — 2 anomaly, which is mismatch between the theoretical and the experimental values 
of the anomalous magnetic moment of muons, provides a sensitive probe of new physics. The 
supersymmetry (SUSY) is one of candidates for new physics models to solve the muon g — 2 
anomaly The minimal supersymmetric extension of the standard model is called the MSSM. 
The MSSM contains superpartner of the SM particles, called superparticles. If masses of the 
superpartners relevant to the muon g — 2 are ^(100 GeV), the SUSY contributions to the muon 
g — 2 become sizable and the anomaly can be solved. 

There are two representative SUSY contributions to the muon g — 2, the chargino and the 
neutralino contributions, respectively. Particularly, when the discrepancy of the muon g — 2 is 
dominated by the neutralino contributions, it is found that the discrepancy can be explained if 
only the Bino and sleptons are light. This setup is "minimal" to solve the muon g — 2 anomaly. 

In this dissertation, we discuss the "minimal" SUSY models where only the Bino and slep¬ 
tons are light, while the other superparticles are decoupled. When masses of relevant super¬ 
particles are ^(100) GeV and the left-right mixing of the sleptons are large, the SUSY contri¬ 
butions to the muon g — 2 are enhanced. In this model, size of left-right mixing is important. If 
it was allowed to be arbitrarily large, the sleptons could be extremely heavy while keeping the 
SUSY contribution, thereby escaping any collider searches. However, it is constrained by the 
vacuum stability condition of the slepton-Higgs potential. Therefore, there are upper bounds 
on masses of sleptons and Bino. 

The searches for the sleptons rely on the slepton mass spectrum. When the slepton soft 
masses are universal, the upper bound on the smuon mass becomes 330 (460) GeV in order 
to solve the muon g — 2 anomaly at the Icr (2cr) level. It is within the reach of LHC and ILC. 
Further, light staus might affect the Higgs coupling to di-photon, which is sensitive to new 
physics. We also study future prospects of the stau which contributes to the Higgs coupling. 
If the stau is heavier than the smuon, the bound can be as large as 1.4 (1.9) TeV. Such non- 
universal slepton mass spectrum generically predicts too large LFV/CPV We show that the 
models are expected to be probed by LHC/ILC and LFV/CPV complementarily in future. 
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Chapter 1 
Introduction 


1.1 Overview 

The Standard Model (SM) of the elementary particle physics is a successful model which 
describes almost all phenomena of the elementary particles. Only the Higgs boson was the 
last piece of the particle contents of SM for a long time. On 4th July 2012, the ATLAS and the 
CMS collaborations claimed that the Higgs like new boson was discovered at last. If this boson 
is the SM Higgs boson, this discovery means that the matter content of SM is completed. 

However, we physicists do not consider the SM as a fundamental theory of elementary 
particles. In fact, new physics models are suggested from some viewpoints, for example, by 
the Hierarchy problem, absent of the dark matter, the dark energy, and the mechanism which 
generated current baryon asymmetry of the Universe, the mismatch between the theoreti¬ 
cal and the experimental values of the anomalous magnetic moment of muons (muon g — 2 
anomaly [5-7]), and so on. We are in the middle of a trip toward the ultimate theory, above 
the suggestions and/or problems are useful guidelines to clarify new physics. 

Particularly, the muon g — 2 anomaly is promising because it provides a sensitive probe 
for new physics contributions. New particles which have masses of ^(100) GeV are suggested 
to explain the discrepancy, and they are within kinematical reach of future collider experi¬ 
ments, such as the Large Hadron Collider (LHC) and the International electron-positron Lin¬ 
ear Collider (ILC). Furthermore, new muon g — 2 experiments are planed at Fermilab [8] and 
J-PARC [9], and the statistical significance is expected to become 7-8 cr in next 3-5 years if 
the central value of experimental would remain the same. In other words, if the muon g — 2 
anomaly is true, it can be confirmed at > 5(j level. Hence, it is important to investigate new 
physics models which are motivated by the muon g — 2 anomaly, just now. In this dissertation, 
we regard the anomaly as basic guideline to search for new physics. 

The supersymmetry (SUSY), which is the symmetry between fermions and bosons, is a 
good candidate for new physics models to solve the muon g — 2 anomally.^ The SM can be ex- 

^ In fact, as we will see in Sec. 3.1, SUSY is often motivated by the hierarchy problem. We consider the muon 
g — 2 anomaly as particularly important problem in this dissertation. Since SUSY can naturally provide a solution 
to the anomaly, we take SUSY models as good candidates for new physics models. For a detail, see Sec. 3.3 
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tended to supersymmetric models, and this model includes new particles, which is associated 
with SUSY, called SUSY particles (superparticles). The minimal supersymmetric extension of 
SM, called the minimal supersymmetric standard model (MSSM), has scalar quarks (squarks) 
and scalar leptons (sleptons) as bosonic partners of the SM fermions. Also the fermionic part¬ 
ners of the Higgs and gauge boson are called the Higgsinos and the gauginos. As we will see 
in Sec. 3.3, if the superparicles couple with the muon, such as smuons, electroweak gauiginos, 
and Higgsinos are around ^(100) GeV, the SUSY contributions to the muon g — 2 become 
sizable and the muon g — 2 anomaly is solved. 

Superparticles have been searched for at the LHC. Since none of them has been discovered, 
colored superparticles, such as gluinos and squarks, are considered to be heavier than about 
1 TeV. Further, the Higgs boson mass of 126 GeV suggests that the stops are as heavy as — 
10)TeV [10-13] and/or have a large trilinear scalar coupling [14]. These two facts contradict 
with the suggestion from the muon g — 2 anomaly, which relevant superparticles should have 
masses of ^(100) GeV. In fact, some of the representative SUSY-breaking mechanism, such as 
the minimal supergravity (mSUGRA) and the gauge mediation (GMSB) cannot simultaneously 
satisfy (i) explaining the 126 GeV Higgs, (ii) avoiding the constraints from LHC experiment, 
and (hi) solving the muon g — 2 anomaly [15]. From the viewpoint of model building. Some 
of extended MSSM have been proposed by some theoretical groups. 

In this dissertation, we take model-independent approach. The above conflict implies that 
the colored superparticles are heavy enough to explaine the Higgs mass and avoid current 
LHC bound, while the superparticles relevant to the muon g — 2 are light (^(100)GeV). The 
candidates of light superparticles depend on dominant SUSY contributions to the muon g — 2. 
There are two representative SUSY contributions as follows, 

(1) Chargino-muon sneutrino contribution (see the left panel of Fig. 3.3), 

(2) Neutralino-smuon contribution (see the right panel of Fig. 3.3). 

When the discrepancy of the muon g — 2 is explained the chargino-muon sneutrino contribu¬ 
tion mainly, the Bino, Wino, Higgsino, and sleptons should be light. On the contrary, if the 
neutralino-smuon contribution are dominant, the Bino and sleptons must be relatively light. 
In order to test whether the SUSY contributions is the origin of the muon g — 2 anomaly, it is 
important to investigate the phenomenogy of these two cases as a first step. 

The case (1) has been studied in Ref. [4] . Authors considered the muon g — 2 motivated 
SUSY models, where only the Bino, Winos, Higgsinos, and sleptons are light, while the other 
superparticles are decoupled. In this case, the searches for (Wino-like) electroweak gaugi¬ 
nos are promising. It has been shown that they can be discovered in near future in most of 
the parameter regions, where the SUSY contributions to the muon g — 2 are dominated by 
chargino-muon sneutrino diagrams. 

The case (2) is "minimal" setup to explain the muon g — 2. In the case (1), the Wino is 
also required to be light, because the chargino contribution to the muon g — 2 is suppressed 
in large Wino mass regions. On the other hand, the neurtalino contribution does not depend 
on the Wino mass, and its phenomenology has not studied sufficiently yet. 


1.1. OVERVIEW 
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Figure 1.1: Mass spectrum of the "minimal" SUSY model, where only the Bino and the sleptons 
are light, while the other superparticles are decoupled. 


In this dissertation, we discuss the phenomenology of the "minimal" SUSY model, in which 
the following superparticles are light, 

B, 4, (1.1) 

while the other superparticles are decoupled,^ as seen in Fig. 1.1. Here, B denotes the Bino 
and i is the sleptons. This study is based on our study [1-3]. The model parameters are 


Ml, m|^. 


m~ 


m~ 


( 1 . 2 ) 


where Mi is the Bino mass, and m~ and m~ are soft SUSY-breaking masses of the left- and 
right-handed sleptons, respectively. m~ is off-diagonal components of the slepton mass ma- 

(■LR 

^ In this dissertation, masses of the heavy superparticles set to 30TeV for simplicity. The results of this 
dissertation are almost independent on their mass spectrum except for the Wino. When the Wino mass are 
< 1 TeV, the chargino-muon sneutrino contribution causes a large theoretical uncertainty. 
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Vacuum meta-stability 
(Sec. 4.1.3) 



U ' B 'in 

-»— *\/V\/X/\/\/\> —•- 

Muon g — 2 
(Sec. 2.2, 3.3) 



Universal 

mg' = mji — rrir 

LHC/ILC 
(Sec. 4.2.3, 4.3) 


Non-universal 
mg' = mji <C mv 

LFV/CPV 
(Sec. 4.2.3) 


Figure 1.2: Overview of this dissertation. 


trices (left-right mixing of the sleptons). Naively if masses of Bino and sleptons are order of 
100 GeV and the left-right mixing is large, the muon g — 2 anomaly can be solved. 

In this model, size of the left-right mixing is important. If it was allowed to be arbitrarily 
large, the sleptons could be extremely heavy while keeping the contribution to the muon g — 2, 
thereby escaping any collider searches. However, too large left-right mixing spoils stability 
of the electroweak vacuum. Therefore, the slepton masses are bounded from above by the 
vacuum meta-stability condition. We will discuss future prospects of slepton searches based 
on the mass bounds. The overview is summarized in Fig. 1.2. 


1.2 Outline of the dissertation 

The outline of this dissertation is as follows. In Chapter 2, we review the electromagnetic 
dipole moments and the Higgs coupling. First, we summarize the electromagnetic form fac¬ 
tors in Sec. 2.1. They induce the magnetic and the electric magnetic moments, respectively. 





1.2. OUTLINE OF THE DISSERTATION 
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The lepton flavor-changing version corresponds to the lepton flavor violating processes. In 
Sec. 2.2, we review the anomalous magnetic moment of the muon, which has discrepancy be¬ 
tween the theoretical and experimental values. Then, the lepton flavor and CP violations are 
highly suppressed in SM. In SUSY models, they are sensitive to off-diagonal components of the 
slepton mass matrices, and are enhanced in case of non-universal slepton masses. Sec. 2.3 and 
2.4 give a review of the lepton flavor violation and the electric dipole moment, respectively. 
The Higgs couplings are also promising. Particularly, since the Higgs coupling to di-photon is 
from loop diagrams, it is sensitive to new physics contributions, such as staus. In Sec. 2.5, we 
review the Higgs coupling to di-photon. 

In Chapter 3, we review the supersymmetry, which is symmetry between fermions and 
boson. In Sec. 3.1, we introduce the supersymmetry. Then, we briefly review the minimal 
supersymmetric standard model (MSSM) in Sec. 3.2. New particles included in MSSM are 
called the SUSY particles (superparticles). If they have the masses of TeV scale, it is expected 
that they affect some observables indirectly. Particularly, observables corresponding to the 
electromagnetic and the weak interactions are sensitive to new physics contributions. In this 
dissertation, we focus on the observables in which the non-colored superparticles are expected 
to have large contributions. The SUSY contributions to the muon g — 2, the Lpy the electron 
EDM, and Higgs coupling to di-photon are reviewed in Sec. 3.3, 3.4, and 3.5, respectively. 
Finally, we discuss the current status for superparticle searches in Sec. 3.6. 

Chapter 4 is the main part and based on Refs. [1-3]. Future prospects for the slepton 
searches are discussed. First, we introduce the "minimal" SUSY model in Sec. 4.1. Then, we 
show that the slepton masses are bounded from above by vacuum mata-stability condition 
in slepton-Higgs potential. The searches for sleptons depend on slepton mass spectrum as 
follows, (1) universal mass spectrum, = m^, and (2) non-universal mass spectrum, 

^ respectively. The former case is discussed in Sec. 4.2.2 and 4.3. The latter 
case is studied in Sec. 4.2.3. 

Chapter 5 is devoted to conclusion. The vacuum decay rate is summarized in Appendix 
A. The corrections to the Bino-muon-smuon coupling are discussed in Appendix B. The cross 
section of slepton productions is summarized in Appendix C. 
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Chapter 2 
Foundation 


This chapter is a review of measurements which are sensitive to new physics contributions. 
In this dissertation, we pick up the anomalous magnetic moment of the muon (muon g — 2), 
the charged lepton flavor violation (cLFV), the electric dipole moment (EDM), and the Higgs 
coupling to di-photon. First, we start from the dipole contributions of leptons, generated by the 
electromagnetic form factor of leptons in Sec. 2.1. They contribute to the muon g—2, cLpy and 
EDM, respectively. In Sec. 2.2, we review the muon g —2, which has a discrepancy between the 
experimental and the theoretical values. The discrepancy is expected to be confirmed by future 
muon g — 2 experiments. Although processes of LEV and CP violation (CPV) are extremely 
suppressed in standard model, some of new physics models might induce large LEV and CPV 
They are also expected to be probed in future experiments. In Sec. 2.3 and 2.4, we review LEV 
and EDM. For a detailed review of the muon g — 2, LFy and EDM, see textbook [16]. 

The Higgs couplings are also promising to probe new physics contributions. Particularly, 
since the Higgs coupling to di-photon is loop-induced, it is sensitive to new physics contribu¬ 
tions. Sec. 2.5 is a review of the Higgs coupling to di-photon. 


2.1 Electromagnetic Dipole Moments 

The dipole contributions of leptons, which induce the electric and magnetic dipole moment 
and flavor-changing processes, not only have been precisely measured, but also are sensitive to 
new physics contributions. In this section, we summarize the contributions of electromagnetic 
dipole moments in field theory approach. 

2.1.1 Electromagnetic Form Factors 

First, we consider the matrix element of the electromagnetic current between 

initial and final states of a spin 1/2 fermion £, with momenta p and p', respectively. Here, we 
define a transfer momentum as q= p' — p. The matrix element is written as 

(£(pO|J^“|((p)) = eQpl^{p')r^u^ip), (2.1) 
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where is the electric charge of the lepton,^ and and are Dirac spinor fields of the 
lepton. The matrix has the general Lorentz structure as [16]^ 

F 

T;, = ~ - 2m^%)r5, (2-2) 

where is the gamma matrices, = (i/2)[y^, y^], and F;(q^) (i = 1,2,3,A), are called the 
charge, anomalous magnetic dipole, electric dipole and anapole form factors,^ respectively. 
The form factors become real since the current is required to be hermitian. 

The static charge and dipole moments are defined at = 0 as follows. 


1^1 (0) = 1 • electric charge, 

1 ^ 2 ( 0 ) = • anomalous magnetic moment, 

d( 

P 3 W = — • electric dipole moment, 
e 


(2.3) 

(2.4) 

(2.5) 


To lowest order, fi(0) = 1 and ^ 2 , 3 ( 0 ) = 0- Non-vanishing contributions to F 23 are generated 
by radiative levels. 

The dipole contributions are also evaluated in effective Lagrangian approach. The effective 
Lagrangian, which provide to ^ 2 , 3 ( 0 ) is given by 


-^eff 


dpAv ~ d^jA^, 




( 2 . 6 ) 

(2.7) 


where F^^ is the field strength and A^ is the electromagnetic field. The interactions in Eq. (2.6) 
are called dimension 5 operators. They are generally not allowed in classical Lagrangian since 
they spoil renormalizability at quantum field theory. Nevertheless, they arise at the radiative 
levels, if no symmetry forbid them. Both the anomalous magnetic and electric dipole moments 
must be finite and are calculable in principle in term of other parameters of the theory. In 
Sec. 2.2 and 2.4, we will summarize the current status and the future prospects for them. 


2.1.2 Transition dipole moments 

Flavor-changing transition amplitudes between distinct fermions are composed of flavor off- 
diagonal part of electromagnetic current and induce —>■ ij + Y processes. They can be 

^ We define e > 0 and Qf T 1 for leptons (anti-leptons). 

^ This definition is slightly different from Ref. [16]. Correspondences of [16] to our definitions are 

E;, ^ eQi ■ r^, 

where the left hand sides are from [16], our definitions are the left hand sides. 

^ The anapole form factor corresponds to parity violations and is generated by electro weak loop physics. 

In this dissertation, we will not discuss anapole induced interactions. 
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parameterized in analogy with Eqs. (2.1) and (2.2) as follows [16] 

(^j(pOI^^”l^i(p)) = (2.8) 

- q/^Pv) f + rs^MoC^")] 

+io-^vq'' • (2-9) 

The matrix is composed of two parts. The first two form factors F^g and F^g correspond to 
chirality-conserving processes at = 0. They are part of dimension 6 operators in effective 
Lagrangian. For example, they mainly contribute to jji —>• eee process. In this dissertation, we 
are interested in scenarios in which the contributions to the dipole moments are dominant, 
thus we will not discuss the dimension 6 operators."^ 

The two form factors in second line of Eq. (2.9), F]]^ and F^^, change the chirality. They 
contribute to flavor-changing processes analog of magnetic and electric dipole moments whose 
contributions are included in dimension 5 operators. They correspond to —>■ + y decays, 

e.g., jJL —» ey. In this dissertation, we assume that dipole contributions are dominant. 

Flavor-changing processes can be also calculated by effective Lagrangian approach. The 
effective Lagrangian which give rise to F^-] and F]^^, is represented as [17] 

-^eff = ^f^o-^.([F^-F^JPi + [F^+F^JpJf^.F'^^+h.c., 
nif - r „ ^ 

= e [AijPi +aIPr) +h.c., (2.10) 

where we define the coefficients andA^^. asA^^. = i.P'ii — Fmi')/^ e,i andA^^. = (f^i 
respectively. Further, . is defined and 0 is assumed for simplicity. 


2.2 Muon g — 2 

The form factor F2(0) in Eq. (2.4) equals to the anomalous magnetic moment of the lepton a^, 
as mentioned in Sec. 2.1.1. Particularly, the muon anomalous magnetic moment, is called 
the muon g — 2. It gives a sensitive probe for new physics contributions as follows [16], 

Experiment: Muons can be copiously produced in a fully polarized state and the lifetimes are 
long enough to precisely measure their precession frequencies in a magnetic field [5]. 

Theory: The SM prediction has been precisely evaluated by several group [18], and the muon 
is heavier than the electron enough to be sensitive to new physics contributions. 

In this section, we briefly review the theoretical prediction and the experimental value of the 
muon g — 2. 

^ Note that the dimension 5 operators also contribute to such as /t ^ eee process even if we neglect the 
dimension 6 operators. 
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2.2.1 Definition 


The g-value, is a dimensionless quantity which characterizes the magnetic moment. It is 
proportional to the spin magnetic moment T^spin, and is defined as 


P 


Spin 


Sf X T— S , 


2 m 


/ 


( 2 . 11 ) 


where m^, eQ^, and S are the mass, electric charge, and spin of the particle /, respectively. 
The quantity jJLg = e/2m^ is called the Bohr magneton. The interaction Hamiltonian 

between the particle and the magnetic field B is expressed as 

Wint=>-7?spin-^- (2.12) 

The g-value for leptons is predicted as 2 in Dirac theory. However, radiative corrections 
shift it slightly. The shift from 2 is known as the anomalous magnetic moment 


=2(1 + a^). 


(2.13) 


In particular, is called the muon g — 2. 

Eq. (2.13) can be checked by field theory based approach in Sec. 2.1. First, we consider a 
lepton scattering from a static vector potential = (0,A(,i). The scattering matrix element ^ 
is evaluated as [19] 


^ = eQi 


Fiiq. )ri + i 




2 m( 


■CTivq 




(2.14) 


where is the Fourier transform of A^,) and the electric dipole form factor F 3 is neglected. 
The spinor u^(p) is expressed as^ 






- / (non-relativistic limit). 




(2.15) 


Then, we extract linear contributions in from Eq. (2.14) as 

J( ^ -(2m,)eQ,r + '^' 2 ( 0 )] ) ^ 0)’ 

where is the Fourier transform of the magnetic field. 


(2.16) 


^ The representation obeys in Ref. [19]. 
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Figure 2.1: (A) The Feynman diagram which contributes to the anomalous magnetic moment 
of the muon. (B) The hadronic vacuum-polarization contributions to the muon g — 2. (C) The 
hadronic light-by-light contributions to the muon g — 2. 


The amplitude ^ can be interpreted as the Born approximation to the scattering of the 
electron from a potential (For detail, see textbook [19]). The interaction Hamiltonian which 
corresponds to such potential is given by 

Hint = -(7?) • 'b, (2.17) 


where 

(TT) = 2 [fi(0) + F2(0)] X (2.18) 

2m^ 2 

The factor /2)^ can be interpreted as the spin of the leptons, S . Comparing Eq. (2.18) 
with Eqs. (2.3) and (2.4), the coefficient 2 [Fi(0) + F2 (0)] becomes 

2[Fi( 0) + F2(0)] =2(l+a,) = g,. (2.19) 

It is just the g-value. 

2.2.2 The Standard Model prediction of the muon g — 2 

The SM prediction of the muon g — 2 has been precisely evaluated through many efforts of 
several groups of theorists. Fig. 2.1 (A) shows the Feynman diagram which contribute to the 
muon g — 2. The theoretical uncertainty reaches to about 0.4 ppm. In this section, we briefly 
review the SM prediction of the muon g — 2. 

QED Contribution 

The quantum electromagnetic dynamics (QED) contributions are the dominant contributions 
to the muon g — 2 (99.993%), and come from the diagrams with leptons and photons. They 
have been calculated analytically up to ^(a^), and recently the full five-loop (^(a^)) contri¬ 
bution has been calculated [20,21]. The present QED value is reported as 

a^(QED) = (11658471.8951 ± 0.0080) X 10“^°, 


( 2 . 20 ) 
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where the uncertainties are from the lepton mass ratios, the ^(a^, a^) terms. The value of a 
is taken from atom, a“^(Rb) = 137.035999049(90). The current uncertainty from QED 
contributions is far below the experimental error and is less important for the discrepancy 
between the theoretical and the experimental value of the muon g — 2. 


Electroweak Contribution 

The electroweak (EW) contributions to the muon g — 2 are from diagrams with the W, Z, and 
Higgs bosons. The one-loop contribution is evaluated as 


^i-ioop(Ew) 

/ i . 


-v/2 871^ 


+ + ^log = 


G, K ”i 


ml 


dx 

0 1 


2x^(2-x) 
X -I- -|x^ 


19.48 X 10 


-10 


( 2 . 21 ) 


where Gp = 1.16637(1) x 10“^ GeV is the Eermi constant, and sin^ = 1 — m^/m^ ~ 0.223 is 
the Weinberg angle. It was calculated by several groups [22-25], after the Glashow-Weinberg- 
Salam model was shown to be renormalizable. In one-loop level, only W and Z boson contri¬ 
butions are relevant at a measurable level. On the other hand, the Higgs boson contribution 
is ^(10“^^), i.e., negligible, due to the small Yukawa coupling of the Higgs boson to muons. 

The EW contribution was calculated at the two-loop level, which is negative [26-29], with 
including leading log effects at three-loop level. Recently, it has been re-evaluated using recent 
LHG result of the Higgs mass [30] . The total EW contributions to the muon g — 2 are 

a^(EW) = (15.36 ± 0.1) X 10“^°, (2.22) 

where the uncertainty comes from hadronic effects in the second-order EW diagrams with 
quark triangle loops, along with unknown three-loop contributions [29,31-33]. The leading 
logs for the next-order term are shown to be small [29,30] . 


Hadronic Contribution 

The hadronic contributions are from QCD interaction. They are composed of the two types 
as follows: (1) Hadronic Vacuum Polarization contribution (HVP), and (2) Hadronic Light-by- 
Light contribution (HLbL), respectively. HVP originates from hadronic quantum corrections 
in the photon propagator. It is shown in Pig. 2.1 (B). The lowest order of HVP is included at 
the higher order correction of HVP and the HLbL are of 
Unlike the QED contribution, HVP cannot be calculated from original QCD, because hadronic 
physics in low energy, i.e., non-perturbative QCD regime, contribute to HVP Nevertheless, 
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thanks to analyticity of the vacuum polarization, the leading effect of HVP can be evaluated 
via the dispersion relation as [16] 


a^(HVP-LO) = 


r 00 


4n^ 


ds K{s) cr‘’(s). 


■^hadrons ’ 


Ami 


(2.23) 


where cr°(s)g+g-^hadrons is cross sections for e^e —> hadrons. They are determined by mea¬ 
sured values as input. The kinematical function K(s) is given by [16] 


K{s) 


X 




--xMogx, 

1 — X 



1 + 


S 


log(l+x)-X + y 


(2.24) 

(2.25) 


In the lowest order, two groups have analyzed by using e"'"e —> hadrons data [7,34] as 

a^(HVP-LO) = (692.3 ± 4.2) X 10“^° [34], (2.26) 

a^(HVP-LO) = (694.91 ±4.27) X 10“^*’ [7]. (2.27) 

The most recent result is evaluation of the next-to-leading order [7], which can be also 
calculated from dispersion relation as 

a^(HVP-HO) = -(9.84 ± 0.07) X 10““. (2.28) 

Another hadronic contribution originates though a light-by-light scattering process. It is 
shown in Fig. 2.1 (C). HLbL cannot be determined from data, unlike HVP It must be calcu¬ 
lated using low energy hadronic models that correctly reproduce property of QCD. In this 
dissertation, we quote two results as follows: 

a^(HLbL) = (10.5 ± 2.6) X 10““, (2.29) 

a^(HLbL) = (11.6 ± 4.0) X 10““. (2.30) 

The former is obtained in Ref. [35], the latter is from Ref. [36,37]. 


2.2.3 Current Status of the muon g — 2 

The complete Standard Model prediction of the muon g — 2 is 


a^(SM) = a^(QED) ± a^(EW) ± a^(HVP ± HLbL). 


(2.31) 
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Contributions 

Value (xlO 

QED 

11658471.8951 ±0.0080 

EW 

15.36 ±0.1 

HVP (LO) [34] 

692.3 ±4.2 

HVP (LO) [7] 

694.91 ±4.27 

HVP (HO) [7] 

-9.84 ±0.07 

HLbL [35] 

10.5 ±2.6 

HLbL [36,37] 

11.6±4.0 

Total SM [34] 

11659180.2 ± 4.2h.lo ± 2.6h_ho ± 0.2„,her(±4.9,oJ 

Total SM [7] 

11659182.8 ± 4.3h.lo ± 2.6h ho ± 0.2„,her(±5.0,„J 

Discrepancy [34] 

28.7 ±8.0 

Discrepancy [7] 

26.1 ±8.0 


Table 2.1: Summary of the Standard Model contributions to the muon g — 2. The lowest- 
order hadronic contribution and the hadronic light-by-light contribution is quoted by using 
two recent evaluations. We take the result in Ref. [7]. 


As for the hadronic contribution, we take the result in Ref. [7] Combining Eqs.(2.20), (2.22), 
(2.27), (2.28), and (2.29), the SM value is 

a^(SM) = (11659182.8 ± 4.9) X 10“^°. (2.32) 

On the other hand, the muon g — 2 has now been measured to a precision of 0.54 ppm by 
the experiment E821 at the Brookhaven AGS [5,6]. The values is summarized as 

a^(exp) = (11659208.9 ± 6.3) X 10”^°, (2.33) 

and the discrepancy between the experimental and the theoretical value becomes 

Aa^ = a^(exp) - a^(SM) = (26.1 ± 8.0) x 10“^°. (2.34) 

It is found that the discrepancy is as large as 3cr level. In this dissertation, we call this problem 
the muon g — 2 anomaly, and assume that the discrepancy suggests new contributions from 
new physics beyond the SM. 

2.2.4 Future Prospects of the muon g — 2 

The theoretical uncertainty on the standard model value is about 0.4 ppm, as seen in Eq. (2.32) . 
We summarize the standard model contributions in Tab. 2.1. Erom this, it is found that the 
theoretical uncertainty by the hadronic contribution is dominant. The lowest-order hadronic 
contribution is determined from the ^ hadrons cross section data. It will be improved 


^ As for HLbL, Ref. [35] is quoted in Ref. [7]. 
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Figure 2.2: The current status and the future prospect for the experimental and the theoreti¬ 
cal values of the muon g — 2. The figure provide comparisons between a^(exp) and a^(SM). 
"DHMZ" is from Ref. [34], "HLMNT" is from Ref. [7]. "SMXX" means the same central value 
with a reduced error as expected by the improvement on the hadronic cross section mea¬ 
surement. "BNL-E821 04 ave." is the current experimental value, and "New (g — 2) exp." is 
the same central value with a reduced error as planned by the future (g — 2) experiments at 
Fermilab [8] and J-PARC [9]. This figure is from Ref. [38]. 


by new data from Novosibirsk and BESIII [38] . The theoretical uncertainty is expected to be 
reduced by a factor 2, i.e., ~0.2 ppm. 

The hadronic light-by-light contribution might be also improved. Two approaches are 
planed to reduce the uncertainty. One is the measurements of y* physics at KLOE [39,40] and 
BESIII. Eor example, the reaction e~^e~ —> e~^e~Y*Y* —> e~^e~T[^ gives a useful information for 
the form factor which is the model parameter of HLbL. The other is the calculation on 

the lattice, which could produce a meaningful result by 2018 [38]. The goal is to compute 
HLbL to 10 % accuracy, or better during the next 3-5 years. ^ 

The new g — 2 experiments are planed at Eermilab [8] and J-PARC [9] . The experimental 
uncertainty is expected to be reduced about 0.1 ppm. Pig. 2.2 shows the current status and the 
future prospects for the theoretical and the experimental values. It is found that the statistical 

^ Improvement for HVP is also planned by the lattice QCD. It is expected to be reduced to 1 or 2 % within the 
next few years. 
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significance would become 7-8 a in near future if the central value would remain the same. 
Thus, if the muon g — 2 anomaly is true, it can be confirmed at more than 5 cr levels. 

The sensitivity to the muon g — 2 anomaly will be improved during the next 3-5 years. 
Therefore, detailed investigation new physics models which are motivated by the muon g — 2 
anomaly are important, just now. Several models have been invented to solve this anomaly. 
The supers)mimetry can be a solution, if the masses of new particles couple to the muon, 
are order of 100 GeV. In this dissertation, we concentrate on the supersymmetric standard 
model. The new contributions from supersymmetric models are given in Sec. 3.3. Searches 
for supersymmetric models which solve the muon g — 2 anomaly are discussed in Chapter. 4. 


2.3 Lepton Flavor Violation 

New physics contributions to flavor changing neutral current (FCNC) processes tend to be 
much larger than the SM predictions. Particularly, muon FCNCs are generically enhanced 
when the muon g — 2 discrepancy is explained by new physics contributions. Therefore, the 
size of FCNC have been strongly constrained from lepton flavor experiments. In this section, 
we summarize theoretical foundations and status for the experimental searches. 


2.3.1 Motivation 


Particle flavors have often provided hints of new physics that was not accessible at the stage. 
For example, the Glashow-Iliopoulos-Maiani mechanism (GIM mechanism) [41] was sug¬ 
gested to explain the absence of quark FCNCs, and predicted the existence of charm quark, 
when only three quarks were discovered yet. Kobayashi and Maskawa hypothesized three 
generations of quarks [42] to explain the CP violation observed in neutral Kaon mixing and 
decays, and after predicted large CP violations in B meson (or extremely heavy quark, i.e., the 
top quark) [43]. This hypothesis was confirmed by discovery of top quark and the B factory 
experiments. Further, the neutrino oscillations, which are evidence tiny neutrino mass, could 
imply the existence of extremely heavy Majorana neutrinos (Seesaw mechanism [44-47]) and 
Leptogenesis [48], a scenario which the CP violation in the neutrino sector might be origin of 
asymmetry between the matter and the anti-matter in the Universe. 

Neutrinos oscillate almost freely among three lepton generations. On the other hand, no 
FCNC processes have ever been observed yet in charged leptons. In the standard model,® they 
are extremely suppressed because neutrino masses are much smaller than the EW scale. For 
example, a process in which a muon converts into an electron and a photon, jU —> ey, can 
proceed through a neutrino oscillation at the one-loop level. The branching ratio of 

11 ^ ey process is evaluated as follows. 


B(iU 


ey) = 


3a 

3 ^ 


S 


m 




MNSJei 2 

77 ? 


® Here, the "standard model" includes neutrino masses. 


(2.35) 
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where Vmns is the Maki-Nakagawa-Sakata (MNS) matrix [49] and m^, are neutrino masses. 
It is found that the branching ratio is negligibly small (< because of tiny neutrino 

masses. In this way, the lepton FCNCs are highly suppressed in SM. 

On the other hand, some of New Physics scenario tends to introduce relatively large 
charged lepton flavor violation (cLFV) since no fundamental symmetry forbid it. Particularly, 
muon FCNCs are more sensitive to new physics contributions. Typical new physics contribu¬ 
tion to the branching ratio of p, —> ey is estimated as 

o I |2 

e 

B{p - er)nv ~ ^2 .V ^ (2.36) 

where A[.jp is a typical mass scale of new physics, and characterizes a size of p-e flavor 
violation. Then, using B(p —> ey) < 5.7 x 10“^^ [50], we can estimate the constraint as 


Anp ^ 600 TeV x 



(2.37) 


If the muon g — 2 anomaly is solved by new physics contributions, masses of new particles 
should be order of 100 GeV. In this case, the size of e^g is constrained as 

e^g < 2 X 10“® for A^p ~ 100GeV. (2.38) 

In this way, muon FGNGs provide a sensitive probe of new physics. Therefore, we concentrate 
cLFV processes of muon in this dissertation. 


2.3.2 Flavor violating muon decays 

In this section, we provide formulae for muon FGNG processes in model independently. We 
assume that dipole moment contributions are dominant. The starting point is the effective 
Lagrangian (2.10) in Sec. 2.1.2. Here, we discuss representative three processes, p —> ey, 
p —» eee, and p — e Gonversion in a Muonic Atom, respectively. 


p —> ey Decay 


One of the most popular LFV processes is the p —> ey decay. The leading contributions proceed 
through the dipole-type operators. Using the effective Lagrangian (2.10), the branching ratio 
of p —> ey is evaluated as [16] 


B(p 


ey) 


r(p —» ey) 487i^a 

r(p ^ all) Gp 




(2.39) 


where A^^ and A^^ ^te Wilson coefficients of the dipole operators. The total decay rate of the 
muon r(p —» all) is well approximated as 


1927i3 


r(p —» all) r(p —> evv) 


(2.40) 
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Process 

Current bound 

Euture sensitivity 

p^ey 

p —»eee 

p — e conversion 

5.7 X 10“^^ (MEG) [50] 
10“^^ (SINDRUM I) [54] 

7 X 10“^^ (SINDRUM II) [56] 

6 X 10-^4 (MEG) [53] 

10“^^ (Mu3e) [55] 

3 X 10"^^ (COMET/Mu2e) [57,58] 

2 xlO"^^ (PRISM/PRIME) [59] 


Table 2.2: Summary of the flavor violating muon decays in Sec. 2.3.3. The left column shows 
muon FCNC processes. The middle column is current experimental limits for each of process. 
The right column represents sensitivities of future experiments. 


jU —» eee Decay 

The dipole operators also contribute to jjl —>■ eee, fJL — e conversion processes via virtual photon 
effects at 7^ 0. In this case, the branching ratio of p, —> eee can be estimated in term of that 
of p —> ey. The ratio between p —> eee and p —> ey processes is evaluated as [16] 


B(p —> eee) a 
B(p —» ey) 371 



11 

T 


0.006. 


(2.41) 


p — e Conversion in a Muonic Atom 

The conversion of a muon captured by a nucleus into an electron has been one of the most 
powerful methods to search for cLFV The ratio of p — e conversion is defined as 


B(pN 


r(piV ^ eN) 
^ r(pN ^ all) ’ 


(2.42) 


where r(piV —> all) is the total decay width. The dipole contributions in this process corre¬ 
spond to that of p —> ey as a function of the mass number (A) and the atomic number (Z), as 
with p —> eee decay. The relation is estimated as [16] 

B(p —> ey) 96n^a 1 428 

-—-—- f2 431 

B(piV^eiV) Gimt 3 X 10^^ BiA, Z) B(A,Z)’ 

r jLi 


where B(A, Z) is a nucleus-dependent factor that includes atomic and nuclear effects, and is 
evaluated as 1.1, 1.8, and 1.25 for Al, Ti, and Pb, respectively [51,52]. 


2.3.3 Current Status and Future Sensitivity 

Signals of cLFV have been searched by many experimental groups. The current sensitivity to 
flavor violating muon decays has reached to about 10“^^. The best upper bound on B(p —> ey) 
is obtained by the MEG experiment at PSI, B(p —> ey) < 5.7 x 10“^^ at 90% CL [50] using 
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~ 1/2 of data taken in summer 2013. The final result of MEG is expected during 2014. 
Further, preparations for an upgrade of MEG experiment (MEG Upgrade) are underway. The 
sensitivity is expected to be » ey) = 6 x 10“^"^ at 90% GL [53] in 3 years of running. 

Searches for jj, —> eee decay are also promising. The current limit is obtained in SINDRUM, 
B(;U —> eee) < 1.0 x 10“^^ at 90% level [54]. Using Eq. (2.41), we find that it is transformed 
to 1.7 X 10“^° in term of » ey decay, and is much weaker than the limit at MEG. New 
measurements is required that an experimental sensitivity is ^(10“^^) at least to be compet¬ 
itive with existing limits and other planned measurements. Recently, the Mu3e experiment 
proposed at PSI. In the first (second) phase of Mu3e, the sensitivity is expected to reach to 
10“^^ (10“^^) [55]. This sensitivity is comparable to that of MEG Upgrade. 

The jJL — e conversion is expected to provide the ultimate sensetivity to cLFV The present 
experimental bound is R{^N —> eiV) < 7 x 10“^^ at 90% GL [56] set by the SINDRUM II 
experiment, in which the muonic gold, i.e., N = Au is used. As with jU —» eee decay, the 
sensitivity to jU — e conversion is expected to be upgraded by two projects. The one is the Mu2e 
experiment at Fermilab, the other is the GOMET experiment at J-PARG. GOMET will follow 
a staged approach, GOMET phase-I and phase-II, respectively. The expected sensitivities of 
COMET phase-I (II) reach to 3 x 10“^^ (3 x 10“^^) at 90% CL with Al [57]. Mu2e experiment 
provides almost same sensitivity, 5.4 x 10“^^ at 90% CL with Al [58]. Further, the PRISM 
project is being developed in Japan. The sensitivity is expected to reach to 2 x 10“^^ at 90% 
CL with Ti [59]. These results are summarized in Tab. 2.2. 

Finally, let us estimate the constraint for mass scale of new physics again, using the ex¬ 
pected sensitivity of PRISM project. Using Eq. (2.43) and B(A.Z) = 1.8 for Al, the sensitivity is 
transformed to 4.8 x 10“^^ in term of jU —> ey decay. In the same way in (2.37), the constraint 
for new physics mass is estimated as 


Anp > 6000 TeV x 



(2.44) 


It is found that sensitivity to new physics scale is improved by order of magnitude in the long¬ 
term future. Similarly, the sensitivity to size of jji-e flavor violation is also estimated as 
e^g < 3 X 10“^° in case that A[.jp = 100 GeV. 

In Sec. 4.2.3, we will discuss the supersymmetric models with non-universal slepton mass 
spectrum, which generically induces large lepton FCNCs. Then, it is found that such models 
are tightly constrained by lepton flavor experiments. 


2.4 Electric Dipole Moment 

There is no experimental evidence for the electric dipole moment (EDM), despite nearly a half- 
century of search. However, laboratories have attempted to find the EDM by new approaches, 
the searches become more important than ever before. In this section, we summarize current 
status of the EDM experiment. Particularly, we focus on the electron EDM, because it provides 
a sensitive probe for new physics models where we will consider in this dissertation. 
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2.4.1 Motivation 

Unless both parity (P) and time reversal (T) invariance are violated, no EDM can exist. The 
interaction Hamiltonian between the lepton I and the electric field E is defined as [16] 

-y, (2.45) 

where S is spin of the lepton. Under a parity transformation, the axial vector S remains 
invariant, but the polar vector E changes sign. On the other hand, under a time reversal 
transformation, E is invariant but S changes sign. In this way, is not invariant under 
P and T transformations. Further, assuming CPT invariance, T violation is equivalent to CP 
violation. Since we know that CP invariance is violated in the decay of neutral Kaon and B 
mesons, the existence of EDM appears quite possible. 

However, EDM is highly suppressed in the standard model. Since the SM has only one CP 
phase in CKM matrix, EDM of the lepton is induced by the CKM phase in the quark sector via 
a diagram with a closed quark loop. Then, a non-vanishing effect appears first at the four-loop 
level [60]. Therefore, the lepton EDM in the SM is extremely suppressed. In electron case, it 
is estimated as 


dCKM < ^Q-38g (2.46) 

This is far from the current experimental sensitivity, ~ cm. 

New physics models genetically contain complex parameters, which induce CP violating 
phase. The new particles which have complex parameters can contribute to EDM at two-loop 
or one-loop levels. It is difficult to justify a suppression mechanism of these CP phases. Hence, 
EDM tends to be enhanced by new physics effects, the observation of a non-vanishing EDM 
would be heralded as a discovery of new physics. 

Let us estimate the new physics contribution. We assume that a typical mass scale of new 
physic is Af^p and parameterizes the size of CP violating phase. Just like the muon FCNCs, 
the typical new physics contribution to the lepton EDM can be naively given by [16] 

m? e 

— tan 0NP, (2-47) 

where the Wilson coefficient of effective Lagrangian is assumed to be ^(1). Then, using 
current bound of the electron EDM < 8.7 x cm at 90% CL [61], we can estimated as 

^NP ~ 240TeV X •\/tan <p^p. (2.48) 

If new particles with ^(100) GeV masses exist, and they could contribute to both the muon 
g — 2 and the electron EDM, the size of CP phase is constrained as 0f^p < 1.7 x 10“^ for 
A^p = 100 GeV. In this way, the electron EDM also provides a sensitive probe for new physics. 
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Material 

Current bound (90% CL, unit of e cm) 

Euture sensitivity 

YbE 

1.05 X 10“^^ [62] 

10“^“ [63] 

Er 


10-29 [ 54 ] 

ThO 

8.7 X 10“^^ [61] 

1 X 10-2®/-y(day) [65,66] 

WN 


10-2° / day [67] 


Table 2.3: Summary of the electron EDM in Sec. 2.4.2. The left column shows atoms or 
molecules which are used to measure the electron EDM. The most severe bound is obtained 
by ThO molecule [61]. The right column is future sensitivities, where "day" represents the 
running time in days. If "day" is ^(10), the future sensitivity will become about 10“^^e cm. 


2.4.2 Current Status and Future Sensitivity 

We cannot directly measure only the contribution of the electron EDM, dg, but search for 
a paramagnetic atom or molecule EDM and then interpret the results in term of the un¬ 
paired electron EDM. Naively, this approaches seem to be impossible due to Schiff's theo¬ 
rem [68] However, this theorem fails when relativistic effects are taken in account as shown 
by Sanders [69,70]. Eurther, for paramagnetic atoms or molecules, since the effective electric 
field is proportional to Z^~^, where Z is atomic number, the EDMs are enhanced. Hence, a 
paramagnetic atom or molecule EDMs with large Z are expected to provide a more sensitive 
probe for new physics effects. 

The bound for electron EDM was obtained by the ytterbium fluoride (YbE) molecules, 
dg < 1.05 X 10“^^e cm at 90% CL [62] until quite recently. The new bound is reported by 
ACME collaboration, and is obtained by the polar molecule thorium monoxide (ThO), dg < 
8.7 X 10“2^e cm at 90% CL [61]. 

The sensitivity is expected to be improved by future experiments. Eor atoms, the future 
sensitivity with the Er atom will reach to dg = 10“^^e cm [64] . Eor molecules, the YbE molecule 
and the WN ion are promising. The sensitivity with the YbE is dg = 10“^° e cm [63]. The 
experiment with the WN ion can probe down to dg = 10“^°e cm, where the systematic limit is 
at the level of 10“^^e cm [67]. These results are summarized in Tab. 2.3. 

Einally let us estimate the constraint for mass scale of new physics again, using the ex¬ 
pected sensitivity by WN ion. In the same way in (2.48), the constraint for new physics mass 
is estimated as 


^NP ~ 7000TeV X (2.49) 

Therefore, it is found that sensitivity to new physics scale is improved by order of magnitude 
in the future, and is comparable to that of muon ECNCs. Similarly, the sensitivity to CP phase 
4>j^p is also estimated as ~ 2 x 10“^° for = 100 GeV. 

^ This is the theorem that the atom or molecule cannot exhibit a linear Stark effect to first order in non- 
relativistic limit. For details, see Ref [ 68 ] or textbook [16]. 
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2.5 Higgs Coupling 

In 2012, a new scalar boson was discovered at the LHC experiment. Its observed properties 
are consistent with the SM Higgs boson within current experimental uncertainty. In addition, 
so far, no extra particles has ever been discovered yet. Therefore, it has been found that the 
SM provides a good description for the elementary particle physics at electroweak scale, in 
both the gauge interaction and the electroweak s)mimetry breaking sector. 

It is important to investigate Higgs couplings to SM particles. If new physics exist near the 
electroweak scale, they affect the couplings indirectly. Particularly, loop-induced Higgs cou¬ 
plings, i.e., the Higgs boson coupling to di-photon, di-gluon or Zy strongly constrain the new 
physics, and they are called the Higgs oblique corrections [71], named after the electroweak 
oblique corrections. In the SM, these couplings are prevented by the gauge symmetry at the 
tree level and are induced by quantum corrections. Therefore, new physics contributions 
might be probed indirectly by measuring the loop-induced Higgs coupling in future. 

In particular, the Higgs coupling to di-photon, which corresponds to the electromagnetic 
interactions, is important. In the SM, it is dominated by one-loop contributions of the elec¬ 
troweak gauge bosons and the top quark. If new physics contain charged particles which 
couple to the Higgs boson, they contribute to the Higgs coupling to di-photon at radiative 
levels. Hence, it is sensitive to the new physics contributions with charged particles. If such 
new particles exist, it is expected that deviations from the SM prediction are observed. 

In this dissertation, we parameterize the Higgs coupling to yy particles and its deviation 
from the SM prediction as 


K 


r 


ghyri^M) 


= 1 + 5k^, 


(2.50) 


where is the Higgs coupling to yy particles and new physics contributions are represented 
by 5k^. The explicit formula for will be given in Sec. 3.5. 

The Higgs coupling to di-photon has been measured at LHC experiment. The mea¬ 
sured uncertainty is 15% (Icr) at ATLAS [72] and 25% at CMS [73], respectively. The re¬ 
sults are consistent with the SM prediction, though they are not yet precise enough to probe 
new physics contributions. In future, the LHC accumulate the luminosity .Sf ~ 300fb“^ at 
^/s = 14TeV, and further upgrade is proposed for .Sf ~ 3000 fb“^ at High-Luminosity LHC 
(HL-LHC). The accuracies are expected to be about 7% (5%) at 300fb“^ (3000fb“^).^° 

Recently, it is argued that the sensitivity to the Higgs coupling to di-photon can be improved 
well, once the international electron-positron linear collider (ILC) will be constructed [74]. 
At LHC, the ratio B(h —>■ —> ZZ*) will be precisely measured, where the Z* means 

Z boson at off-shell. On the other hand, at ILC, the Higgs couplings to Z boson and di¬ 
photon can be measured at (sub) percent levels [75]. The joint analysis of HL-LHC and 
ILC enable us to realize the accuracy of of about 2% [74]. Here, it is assumed that the 
uncertainty of B(h —> YY)/B{h —> ZZ*) is 3.6% from HL-LHC and ILC runs at -v/s = 250GeV 

They are dominated by systematic uncertainties.The accuracies could be improved by reducing them. 
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Figure 2.3: The accuracy of the Higgs coupling to di-photon jc^. The horizontal axis is lumi¬ 
nosity of ILC and the vertical axis is the accuracy of CMS-1 (2) are the future prospects of 
at CMS for .Sf = 3000 fb“^. CMS has presented a less conservative set of estimates, using 
two scenarios. ILC (red) is the accuracy of at ILC and ILC -I- LHC BR ratio (aqua blue) is 
combined with the measurement of B(h —> yy)/B(h —> ZZ*) by ATLAS. Using joint analysis of 
HL-LHC and ILC, the accuracy of will be improved well. This figure is from [74]. 



CMS 250 500 500up 1000 1 OOOup 


and .Sf = 250 fb“^. The direct measurement of at ILC is not so precise that of LHC, due to 
lack of luminosity. If ILC accumulate more luminosity, e.g., if' = 2500 fb“^ at ^/s = ITeV, the 
accuracy of can become 1.9% [74,75] at ILC. the joint analysis can be better than 1% [74] . 
They are so precise that the new charged particles could be sufficiently probed by measuring 
Ky. These discussions are summarized in Table. 2.3. 

In Sec. 4.3, we consider a situation that an excess of is measured in HL-LHC and ILC, 
and consider new physics models that explain the anomalous excess. Although many models 
have been proposed, we focus on the supersymmetry (SUSY), which includes a scalar partner 
of the tau lepton (stau). Then we will study the properties of staus that are responsible for the 
fc^ excess, and discuss whether the stau contribution to can be probed at ILC. 













36 


CHAPTER 2. FOUNDATION 


Chapter 3 

Supersymmetric Standard Model 


The chapter is a review of the supersymmetry. The supersymmetric extension of the stan¬ 
dard model provides a beautiful solution to the hierarchy problem of the standard model. In 
Sec. 3.1, we introduce the supersymmetry motivated by the hierarchy problem. The minimal 
extension model is called the minimal supersymmetric standard model (MSSM). In Sec. 3.2, 
we provide MSSM Lagrangian. The supersymmetric model includes superpartners of standard 
model particles associated with supersymmetry, which is called superparticles. If superpar¬ 
ticles have masses of ^(100) GeV, they are not only discovered by collider experiments, but 
also affect low energy measurements through radiative corrections. Particularly, searches for 
superparticles without color charge are promising in near future experiments. In Sec. 3.3-3.5, 
we summarize measurements which non-colored superparticles could provide large contribu¬ 
tions. Sec. 3.6 is a summary on current status of supersymmetry. The prospects for searches 
of sleptons will be discuss in Chapter. 4. 


3.1 Motivation 

The Standard Model (SM) provides a very successful description of phenomena of the ele¬ 
mentary particles at the electroweak scale rriy^. Although no additional structures have ever 
discovered yet, physicists do not consider the SM as the ultimate particle theory, but a low 
energy effective theory of a fundamental theory, whose typical energy scale is characterized 
by Muy. For example, the reduced Planck scale Mp = — 2-4 x lO^^GeV, which 

characterizes a quantum gravity, is one of important fundamental scales. Once we take above 
standpoint, the Higgs sector includes an "unnaturalness" in the SM. 

This unnaturalness comes from the squared Higgs mass parameter m^. First, we introduce 
the SM Higgs field H with a potential 

V = ml\H\^ + ?i\H\\ (3.1) 

where A is the quartic Higgs self-coupling. The Higgs field has a vacuum expectation value 
(VEV), (H) = / —m^/2A at the minimum of the potential, if A > 0 and < 0. Since the 
VEV is known as about 174GeV, the mass parameter must be order of —(lOOGeV)^. 
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H 



H 


S 

0 

H -^- H 


(a) (b) 

Figure 3.1: The radiative correction to the Higgs squared mass parameter where (a) is 
from a Dirac fermion f, and (b) is from a scalar S. 


The squared Higgs mass parameter receives quantum corrections from diagrams in Fig. 3.1. 
When the Higgs field couples to a Dirac fermion f, the correction as Fig. 3.1 (a) become 

l^fT 

Am^ = + ^(log(M^)). (3.2) 

On the other hand, 

A 

Am^ = + ^(log(M^)), (3.3) 

is from a scalar S in Fig. 3.1 (b). 

In the SM, top quark has strong coupling with Higgs field, ~ 1. However, scalar fields 
other than Higgs does not exist. Hense, we need a fine tuned cancellation between the bare 
mass and the corrections which are proportional to as 

m^(physical) = m^(bare) + Am^, (3.4) 

to realize = —(^(lOO)GeV)^. If Muy = Mp, size of the quantum correction is some 30 
orders of magnitude and is extremely larger than This unnaturalness is known as the 
hierarchy problem [76-80], and is one of motivations introducing models beyond the SM. 

Supersymmetry (SUSY) [81], a symmetry between fermions and bosons, provides a beau¬ 
tiful solution for the problem, and introduces two scalar boson for one Dirac fermion, and 
vice versa. Further, SUSY guarantees that couplings of scalar partners with the Higgs field are 
same as those of Dirac fermions, . Thanks to new contributions of superparticles, 

the quadratic divergence of the Higgs mass can cancel as 

IA P A 

^ ^ ^ ^ ^(log(M^)) - ^(log(M^)). (3.5) 

In this case, the tuning is not only so fine, but also we need not care quantum corrections at 
any energy scales. In this way, the supersymmetric extended models of the SM can solve the 
problem of the SM. Hereafter, we adopt SUSY models as candidates of new physics. 
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Matter and Higgs fields 
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-2/3 

-1/3 
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Dt 
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1/3 

-1/3 

0 

d; 

dl 

Li 

1 

2 

1/2 

0 

1 



Ei 

1 

1 

1 

0 

-1 

7* 


Hu 

1 

2 

1/2 

0 

0 



Hd 

1 

2 

- 1/2 

0 

0 



Gauge fields 


SU(3)g 

SU(2) 

L U(l)y 

B 

L 

spin 1/2 

spin 1 

G 

8 

1 

0 

0 

0 

g 

g 

W 

1 

3 

0 

0 

0 

W 

W 

B 

1 

1 

0 

0 

0 

B 

B 


Table 3.1: Particle contents of the MSSM. 


3.2 The MSSM 

The minimal supersymmetric standard model (MSSM) is the minimal supersymmetric exten¬ 
sion of the standard model. The gauge symmetries of the MSSM are SU(3)cXSU(2)iXU(l)y, 
just like the SM. The MSSM provides rich phenomenology since it includes many new parti¬ 
cles coupled with SM particles and new parameters. In this section, we introduce the MSSM 
superpotential and the soft supersymmetry-breaking terms. 

The superpotential for the MSSM is given by [82] 

Vk = H,L,Ej, (3.6) 

where H^, Q, L, U, D are chiral superfields in Tab. 3.1. Unlike the SM, we need two 
Higgs fields and because of holomorphy of superpotential and anomaly cancellation 
condition. Parameters ^ g are the Yukawa coupling matrices whose indices correspond to 
the generation of the fermions. All of the gauge (color and weak isospin) indices in Eq. (3.6) 
are suppressed. The first term of right hand side, is called jj, term, where 

is the invariant tensor for SU(2)i and indices a,p = 1,2. Similarly, the term H^Q^Uj is 
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Figure 3.2: Feynman diagram for proton decay occurred by the interactions UDD and LQD. 
These interactions cause dangerous contributions to proton decay To forbid such a interaction, 
we have to impose the R parity 


written out as ^uaQ-ipaUj where a = 1,2,3 is a color index which is lowered (raised) 

in the 3 (3) representation of SU(3)c- 

The superpotentilal (3.6) is minimal in the sense that it is sufficient to produce a phe¬ 
nomenologically viable model. The most general gauge-invariant and renormalizable super¬ 
potential includes not only Eq. (3.6) but also the terms [82] 



(3.7) 

y^AB=i = ^X"^^%DjD„ 

(3.8) 


where i = 1,2,3 correspond to generation of the fermion, and B, L are baryon and total lepton 
numbers, respectively. The chiral superfield Q; has B — +1/3, Ui,D^ have B = —1/3 and B — 0 
for all others. The lepton number assignments are L — +1 for Lj, L = — 1 for £; and L — 0 
for all others. Hense, the terms in Eq. (3.7) violate lepton number by 1 unit and those in Eq. 
(3.8) violate baryon number by 1 unit. 

The B- and L- violating processes have not been discovered yet. The most severe constraint 
comes from no observation of proton decay. If the couplings A^, X" were not suppressed, the 
lifetime of proton would becomes extremely short. Eor example, Eig. 3.2 is the Eeynman 
diagram with squark exchange for p~^ —> e~^Ti^ decay. The decay rate is roughly estimated as 




m 


5 

proton 


Zl 

i=2,3 


/Hi A //Hi 





2i.2,3 

Aj 
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1 1 

f 1 TeV^ 
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(3.9) 
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The proton decay is strongly constrained for the Super-Kamiokande experiment [83]. When 
Tp+^g+^o is compared with the experimental constraint, one finds < 10“^^. 

To solve above unnaturalness, we usually impose conservation of the H-parity [84], which 
is a discrete Z 2 symmetry defined as 



(3.10) 


where s is the spin of the particle. If we impose the H-parity conservation in the MSSM, the 
terms in Eqs. (3.7), (3.8) are forbidden. 

The H-parity takes values of = ±1. The P-parity odd particles are called "superparticles". 
If P-parity is exactly conserved, there are no mixing between superparticles (P^ = —1) and SM 
particles (P^ = +1). Further, every interaction in the MSSM contains the even number of 
superparticles. These have three important phenomenological consequences [82]: 

• The lightest particles with P^ = — 1 is called "lightest supersymmetric particle (LSP)", 
and must be stable. If LSP is electrically neutral and interacts weakly with ordinary 
particles, it can be candidate for non-baryonic dark matter that seems to be required by 
cosmological consequences. 

• Each superparticle other than LSP decays into states with the odd number of LSPs. 

• Suparparticles can be produced in even numbers at collider experiments. 

In this dissertation, we concentrate on the MSSM with P-parity. 

We can construct the full supersymmetric Lagrangian using kinetic terms and the super¬ 
potential given in (3.6). However, it is not sufficient to explain phenomena of the Nature. 
The most important constraint is that none of superpartners of the SM particles has been dis¬ 
covered yet. For example, if SUSY is unbroken, the selections must have same mass of 
the electron, 0.511 MeV. A similar discussions can apply to each of the other sleptons, 
squarks, and gauginos, respectively. Such superparticles should be already detected if they 
exist. Hense, SUSY must be a broken symmetry in our vacuum. 

A hint to describe the nature of SUSY breaking, is the hierarchy problem as introduced in 
Sec. 3.1. The hierarchy problem is solved in supersymmetric models because SUSY guarantees 
cancellation of the quadratic divergence to the Higgs mass. If broken SUSY theories still 
provide a solution to the hierarchy problem, the relationship between dimensionless couplings 
that hold in an unbroken SUSY theory must be maintained. Therefore, we are led to consider 
"soft" supersymmetry breaking. This means the effective Lagrangian of the MSSM is written 
in the form [82] 


(3.11) 



where .^svsy is the SUSY invariant Lagrangian and contains the SUSY breaking interac¬ 
tions with coupling parameter whose mass dimension is "positive". In this case, the cancella¬ 
tion of the quadratic divergence is maintained. 
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Then, we provide the soft SUSY breaking terms. They are composed of the mass term of 
the gauginos, the sfermions and the Higgs bosons, and the trilinear term of the scalar bosons 
as follows [82]^ 




-- (MiBB + M 2 WW + Mggg + h.c.) 

+ MijHuQiUj - - (aJijHdL^ej + h.c. 

(m|),. QlQj + Ul] y Lh + (^o) ,j + [<] y + H] y 

+ h.c.)] , (3.12) 


where M^, m~ (/ = Q,L,U,D,E'), and mf, are the masses of the gauginos and the scalar 

f “u,d 

bosons (Higgs bosons and sfermions), ^ ^ are trilinear scalar couplings called A-term, and 
b is quadratic coupling of the Higgs boson called-B term. Note that ^ is not the chiral 
superfields but the scalar fields. Each of ^ ^ is a complex 3x3 matrix. They are in one- 
to-one correspondence with the Yukawa couplings of the superpotential. Each of m~ is a 
Hermitian 3x3 matrix. In this dissertation, we parametrize the A term and B term as 



The soft SUSY breaking terms S^soft introduce many new parameters, which do not exist 
in the SM. There are 105 parameters of masses, phases, and mixings in the MSSM [85] that 
cannot rotate away by redefining the phases and the bases for the fermion supermultiplets. 
These parameters introduce new flavor/CP violations, which might easily contradict with cur¬ 
rent experiments. These problems are called the SUSY flavor/CP problems. Thus, in principle, 
SUSY breaking induces a tremendous arbitrariness in the Lagrangian. 

The size of the flavor/CP violations depends on the origin of the SUSY breaking. In this 
dissertation, we do not enter the detailed models of SUSY breaking, but consider the MSSM 
with the soft SUSY breaking terms as a effective theory at TeV scale. Then, we investigate the 
constraints of low energy measurements to probe the possible type of models. In subsequent 
chapter, we discuss the SUSY contributions to several measurements and introduce the current 
status of SUSY models. 


3.3 Muon g — 2 in the MSSM 

As mentioned in Sec. 2.2, the muon g — 2 has more than 3cj discrepancy between the exper¬ 
imental and the theoretical values. The size of the discrepancy is order of 10“^. It is difficult 

^ In general SUSY models, includes the tadpole terms, such as t'(/>;. They do not occur because there is 
no singlet field in the MSSM. 
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for typical new physics models to explain the discrepancy. First, let us introduce this difficulty. 
Naive new physics contributions to the muon g — 2 are expressed as 



(3.14) 


where a^p and m^p are typical coupling and mass of new particle. 

Since the size of the muon g — 2 discrepancy is comparable with the SM electroweak con¬ 
tributions (2.22), typical masses of the new particles are required to be < if the coupling 
is about a 2 = g^lA-n, where g is the 811(2)^ gauge coupling. However, no such particles have 
ever discovered yet. On the other hand, when typical new particles masses are large enough to 
avoid current mass bound, the new contributions is naively too small to explain the muon g —2 
discrepancy. Therefore, new physics models which can explain the muon g — 2 discrepancy 
should be categorized into the two type as follows: 

(1) The models in which the new particles have small coupling with the muon (^(10“^)) 
and small mass (^(l)MeV), for example. Dark Photon model [86,87]. 

(2) The models in which the new particles have large coupling with the muon (^(1)) and 
large mass (^(100) GeV), for example. Supersymmetry. 

As we will see in Sec. 3.3.3, SUSY contributions to the muon g — 2 are proportional to 
tan 13, which is the ratio of the Higgs VEVs, tan 13 = (H„) / {H^). If tan ^ is order of 10, ^(10“^) 
contributions to the muon g — 2 are realized by masses of ^(100) GeV smuons, which is 
within kinematical reach of future experiments. In this section, we summarize the MSSM 
contributions to the muon g — 2. 

3.3.1 Relevant interactions 

In general, particles coupled with muon affect the muon g—2. In the MSSM, new contributions 
are induced by the superpartners of muon, muon neutrino, SU(2)iXU(l)y gauge bosons, 
and Higgs bosons, respectively. The corresponding interactions are those of gaugino-slepton- 
leptons. In this section, we provide general relevant interactions including case of lepton 
flavor and CP violations. 

First, let us discuss the sleptons. We denote by and the superpartners of leptons 
and neutrinos v^, respectively. In this dissertation, we assume that the slepton mass matrices 
are universal among the flavors, such as in the gauge mediated SUSY-breaking models [88], 
while the Yukawa matrices are non-diagonal in the model basis. After electroweak symmetry 
breaking, the mass matrix of the charged sleptons are written by 



(3.15) 
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where the matrices ^ are written as 



= diag(m~ , m~ , m~ ). 

eiP Mll’ 

(3.16) 

^Ir 


(3.17) 

^Ir 


(3.18) 


Here, A^ is A-term for the sleptons, jtt is the Higgsino mass parameter, and is the lepton 
Yukawa matrix, which is diagonalized by two unitary matrices, U^, L/r as 

U^Y,Ul = diag(7/-^Y^"-^7f^sj. (3.19) 

Here, U^, L/r generically include mixings among the flavors and complex phases, which cause 
LFV and CPV The explicit forms will be given in Sec. 4.2.3. The squared masses of the left- 
and right-handed sleptons m~ and m~ (£ = e,ju, t) are 

^LL ^RR 

m- = m~ 3-ml + Dj (3.20) 

where, Dj^^ are D-terms, Dj = m^ cos2^(r| — sin^ 9^^). 

We diagonalize the mass matrix of the sleptons by 6 x 6 unitary matrix as 

Iji i^LL ^LR I jjt\ — (diagonal). (3.21) 

\Kr Kr’J 


The mass eigenvalues are denoted by (A = 1, • • • 6) and the mass eigenstates are written 
= (i = 1,2,3). 

The mass matrix of sneutrinos is diagonal in the model basis. Further, there are no right- 
handed sneutrinos in the MSSM. Hence, in this dissertation, no sneutrino mixings between 
distinct flavors appear, except for Yukawa interactions. 

Then, we consider the charginos. The mass matrix of the charginos is written as 




M2 

\/2m^y cos P 


V2myY sin ^ 
P 



(3.22) 


It is diagonalized by two unitary matrices, O^, Oj^. We denote the eigenvalues by (A = 
1,2). The mass eigenstates are called the charginos and are defined as 


Zi = (Ol) 




Then = Xal'^ Xar(-^ = 2) forms a Dirac fermion. 


(3.23) 
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Finally, we consider the neutralinos. The mass term of the neutralinos is written as 


where 







uL 




VKJ 


+ h.c., 


(3.24) 




f 

0 

—sin 6y/ cos ^ 
cos cos /3 


0 

M 2 

sin sin /3 
—cos 6y/ sin /3 


—sin cos jS 
rriz cos 6^ cos /3 
0 


sin sin P ^ 
—cosS^sin^ 


(3.25) 


The mass matrix is diagonalized by unitary matrix O^. The mass eigenvalues are denoted 
by myo (A = 1,4) and are called the neutralinos. They are defined as 

AA 


r,= (o^) 


f 


(3.26) 


Then, X^ (^ = 1, • • • , 4) forms a Majorana fermion. 

We are ready to write the relevant interactions. The interaction Lagrangian for gaugino- 
slepton-lepton is written as 


.Sf- 

—^in 


I 


R(e) 


iAX 


P,+N 


L(e), 

iAX ^ 


) Xa^x +^i [cSx^Pr + Xa % + h-c., 


(3.27) 


where the coefficients are 


C’ = -s(0h,,5„. 

cS’ = 

= ■^((0|v)Ai'an®»' + (0|v)M 

= -'/2«tane„(o;)^,[;{y3- 




(3.28) 

(3.29) 

(3.30) 

(3.31) 


3.3.2 Formulae in mass eigenstates 

At leading order, the new contributions to the muon g — 2 are induced by the chargino-muon 
sneutrino and the neutralino-smuon loop diagrams in Fig. 3.3. Each diagram is called the 
chargino-muon sneutrino and the neutralino-smuon contribution, respectively. Computing 
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Figure 3.3: The diagrams of the SUSY contributions to the muon g — 2 in mass eigenstates. 
The left panel of the figure is called the chargino-muon sneutrino contribution (or chargino 
contribution), the right is the neutralino-smuon contribution (or neutralino contribution). 


the Feynman diagrams by ordinary perturbation theory, we can obtain the SUSY contributions 
to the muon g — 2 at one-loop leveP as [89,90] 


Aa^ 


Aa^ 



I^C 


L(e)* 

2AX 







(3.32) 


(3.33) 


where the loop functions. 


Ff(x) = 

" 1 

^2 -F 2x — 6x^ + x^ -F 6x logxJ , 

(3.34) 

(l-x)4 1 

= 

3 

■ [-3 -F4x — x^ - 21ogxJ , 

(3.35) 

2 (1-x)3 

ff(x) = 


— 6x -F 3x^ + 2x^ — 6x^logxj , 

(3.36) 

(i-^r 1 

Ff(x) = 

" 1 

— x^ -F 2x logxj . 

(3.37) 

(1 — x)^ 1 
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Figure 3.4: The diagrams of the SUSY contributions to the muon g — 2 in gauge eigenstates, 
where the diagrams (a)-(e) corresponds to Eqs. (3.38)-(3.42), respectively. The diagram 
(a) comes from the chargino-muon sneutrino diagram, the diagrams (b)-(e) are from the 
neutralino-smuon diagram. 


3.3.3 Formulae in gauge eigenstates 

Even though the formulae in mass eigenstates (3.32) and (3.33) are useful for numerical cal¬ 
culations, they are not suitable for understanding their dependence of the MSSM parameters. 
The parameter dependences are hidden by the electroweak symmetry breaking which causes 
complicated mixings. On the other hand, in weak eigenstates, the structure of one loop con¬ 
tributions become more transparent. The expression in weak eigenstates are equivalent to the 
expansion of where represents SUSY-breaking masses and the Higgsino mass 

jjL. At the leading order of rnyy/nisaft and tan p, they are evaluated as [4,90] 

__ _ (7 f M^ u?' A 

Aa^(W,H,v,) = ^^tan^-/J^,-^^ , (3.38) 

'^471 M 2 E- V. nt£ mij 

^ ^ ~ ay rnlMy^ 

Aai^L,ER,B) = - -^ 

^ An mi mi 

Mi Mr 

Aa^(W,H,pEJ = -^-^tanl3-fj^,-^], (3.40) 

8n M 2 IJL V Mi MiV 


mt 


Mi 


mt 


Mr 


Mf M: 


1 


(3.39) 


^ Strictly speaking, these results include the flavor/CP violating effects which are induced by off-diagonal 
elements of the slepton mass matrix. Since they are subdominant contributions, we neglect them in this section. 
These effects are discussed in Sec. 3.4 and 4.2.3. 
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Mass spectrum 

(a) (b) (c) (d) (e) 

Note to solve the anomaly 

Degenerated case 

y 


Large LR mixing 

y 

Both are light 

W, jU^-decoupled 

y 


W, jU^-decoupled 

y 

IX <0 for a positive contribution 


Table 3.2: Summary on the SUSY contributions to the muon g — 2. The items (a)-(e) in the 
middle column correspond to the diagrams (a)-(e) in Fig. 3.4. The left column shows several 
mass spectrums. The check marks in the middle column denote the dominant contribution to 
the muon g — 2 in the case of left column. 




av 


m. 


8n MjjU 


tan p ■ fi 


N 


av 


m. 


471 M^jJL 


tan 13 ■ fi 


N 



(3.41) 

(3.42) 


where is the muon mass, while ay and a 2 are the fine structure constants of the U(l)y 
and the SU(2)i gauge symmetries, respectively. They can be obtained from the Feynman dia¬ 
grams shown in Fig. 3.4. The arguments in the left hand side of Eqs. (3.38)-(3.42) represent 
the superparticles which propagate in each loop diagram. If one of them decouples, the cor¬ 
responding contribution is suppressed. Eq. (3.38) comes from the chargino-muon sneutrino 
diagrams, and Eqs. (3.39)-(3.42) are from the neutralino-smuon diagrams. 

The loop functions are defined as^ 


/iv(^,y) 




^y 


5 -3(x + y)-Fxy 

(X - l)2(y - 1)2 

-3 + x + y+ xy 


21ogx 


+ 


(x-y)(x-l)3 
2x logx 


+ 


2 logy 


(^-y)(y-i)^ 

2ylogy 


(X - l)2(y - 1)2 (X - y)(x - 1)3 (X - y)(y - 1)^ J 


(3.43) 

(3.44) 


They are satisfy 0 < /c,N(^,y) < 1 and are monochromatically increasing for x > 0 and 
y > 0. When the superparticles have degenerate masses, they become /^(l,!) = 1/2 and 

/jv(l,l) = l/6. 


3.3.4 Structure of SUSY contributions 

As discussed in Sec. 2.2, the muon g — 2 has the discrepancy at more than 3cj levels between 
the experimental and the theoretical values. In order to explain this discrepancy, the new 


^ The functions, and /jv, are reduced from the functions, J 5 and U, in Ref. [90]. 
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contributions to the muon g — 2 are required to be and positive. In this section, 

we discuss the results obtained in Eqs. (3.38)-(3.42) from above viewpoints. The results are 
summarized in Table. 3.2. 

Numerically, the SUSY contributions are evaluated as [4] 


Aa^(W,H,v^) ~ 15x10 

^ 1.5x10“^ 
Aa^(VY,H,^J ~ -2.5x10- 
~ 0.76x10 


f tan f (100GeV)2^ f fc 


V 10 7 M2^ J yi/2j’ 

tan 13 \ f (100 GeV)^ \ f 


10 


nd ml iM.ix \ 

iJ-L Mr' le" / V / y 


N 


g I tan^^ f (100 GeV)^^ f 


{^) i 


M2B 


i/ej’ 


7 tanP'\ f (100 GeV)^^ f ff^ 
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(3.45) 

(3.46) 

(3.47) 

(3.48) 

(3.49) 


where Eqs. (3.38)-(3.42) correspond to Eqs. (3.45)-(3.49), respectively. As mentioned in Sec. 
3.3.3, the arguments in the left hand side of Eqs. (3.45)-(3.49) represent the superparticles 
which propagate in each loop diagram. They are enhanced when tan ^ is large and 
is small. If relevant superparticles which propagate in each loop diagram have masses of 
^(100)GeV and tan/3 = ^(10), they become ^(10“^), which can explain Eq.(2.34). 

Let us investigate parameter dependences of the SUSY contributions in detail. There are 
two representative cases when the discrepancy of the muon g — 2 is explained. The first case 
is the one that the Wino-muon sneutrino contribution (3.38) dominates the SUSY contribu¬ 
tions. This situation is realized when relevant superparticles have nearly degenerate mass 
spectrum, as can be seen in Eqs. (3.45)-(3.49). In fact, the Wino-muon sneutrino contribu¬ 
tion is dominant in wide SUSY models which explain the discrepancy of the muon g — 2 with 
non-decoupling Higgsinos. 

The second case is the one that the pure-Bino contribution (3.39) is dominant. As jji 
increases, it becomes relevant because it is proportional to jU including in left-right mixing of 
the smuon mass matrix. The other contributions are suppressed because they are proportional 
to 1 /jU, which is from the Higgsino propagator. In this case, the muon g — 2 discrepancy is 
explained when both and are required to be ^(100) GeV and tan p = ^(10). 

The other extreme cases also exist. Here, we comment on two cases as follows, (W, Jljd 
and {W, Jli) decoupled cases, respectively. Eirst, when Wino and right-handed smuon are 
decoupled, the Bino-Higgsino-left-handed smuon contribution (3.41) is relatively relevant, 
while the others are suppressed. On the other hand, the Bino-Higgsino-right-handed smuon 
contribution (3.42), which gives positive contribution in p, < 0 case, is relatively relevant, 
when the Wino and the left-handed smuon are decoupled. In these cases, since light smuon 
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and large tan are required to explain the discrepancy of the muon g — 2, the parameter 
regions are extremely limited. In this dissertation, although these cases are also interesting, 
we will not discuss them anymore. 

The Wino-Higgsino-left-handed smuon contribution (3.40) has the same parameter de¬ 
pendence as the Wino-muon sneutrino one, M 2 , ijl, and respectively. As can be seen in 

Eq. (3.45) and (3.47), it can’t be dominate in any parameter region. 

Since the SUSY contributions are proportional to contributions (3.38)-(3.42) 

are rapidly suppressed, as both masses of the smuons increase. In this case, no solution exists 
to explain the discrepancy of the muon g — 2. ^ 

So far, we discussed the SUSY contributions to the muon g — 2 at one loop levels. Some 
of higher order corrections are important in numerical calculations. In this dissertation, the 
corrections that can be as large as or larger than ^(10)% are included and will be discussed 
in Sec. 4.1.2. 

3.4 LFV and EDM in the MSSM 

The MSSM, which explains the muon g — 2 discrepancy, contains light superpartilces coupled 
with the muon. In this case. Lepton Flavor Violation (LFV) and CP Violation (CPV) which are 
sensitive to off-diagonal components of the slepton mass matrices, might be induced. They 
are suppressed if the slepton mass matrices are universal among the flavors, such as in the 
gauge mediated SUSY-breaking models. However, even when the mass matrices are diagonal 
in the model basis, sizable FCNC and CPV are generically induced as long as the diagonal 
components are not equal to each others. In fact, a lot of models have been proposed to 
explain the SM Yukawa couplings, and many of them predict non-diagonal Yukawa matrices 
in the model basis. Even if the sfermion mass matrices are diagonal in this basis, off-diagonal 
components are generated in the fermion mass eigenstate basis, which are obtained by rotating 
the mass matrices of the model basis with unitary matrices. In other words, the super GIM 
mechanism does not work generically unless the sfermion mass matrices are universal. 

In this dissertation, we focus on muon FCNCs and electron EDM. This section is summa¬ 
rized these formulae at leading order. We do not consider the higher order corrections, i.g., 
Barr-Zee type diagrams, because the contributions at one-loop level are not suppressed and 
the heavy Higgs bosons are decoupled in our setup. 

3.4.1 Formulae in r(lj £“7) 

We start to the effective Langangian (2.10). The diagrams which contribute to the ij —»iTy 
process have the same stricture as that of the muon g — 2, as shown in Fig. 3.5. Calculating 

rrijx^^ ~ 100 GeV and tan^ = 40-50 are required. 

^ the muon sneutrino mass is related to mn by the SU(2), symmetry as m? = -I- cos 2/3. 

® Strictly speaking, when the left-right mixing of the smuons is extremely large, ^(l)TeV masses of the smuons 
is allowed. We will discuss the case in Chapter. 4. 
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Figure 3.5: The diagrams of Ij decay. The left panel of the figure is the chargino- 

sneutrino contribution, the right is the neutralino-slepton contribution. 


the Feynman diagrams, the decay rate for —> £. y is evaluated as 

r«7 - <-r) = (|a‘ +<|") . o.so) 

where are the amplitude which chargino-sneutrino and neutralino-slepton contribute to, 
and are given as 
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, (3.51) 

, (3.52) 


and Ap = A^ The loop functions E^^ were given in Eqs. (3.34)-(3.37). 

When j = 2 and i = 1 (i.e., jJL ey decay), dominant decay mode of muons is yT —> 
e“v^Vg. The branching ratio of ;Lt —> ey decay can be calculated in Eq. (2.39). 

Once we can calculate Eqs. (3.51), (3.52), y —> eee and y—e conversion processes are eval¬ 
uated in Eqs. (2.41) and (2.43), if the dipole contributions are dominant. In the MSSM, the 
dipole contributions are proportional to tan p, while non-dipole contributions, which are from 
the dimension 6 operators, are not so. Since large tan p is required to explain the discrepancy 
of the muon g — 2, the dipole contributions are dominant in our assumption. 

Let us comment on higher order contributions, e.g., the Higgs-mediated ECNC [91]. In 
supersymmetric models, the superpotential must be holomorphic. Even though, one can diag¬ 
onalize the mass matrices of fermion at tree level, off-diagonal components of Yukawa matrices 
are generated by radiative corrections, which pick up SUSY breaking effects as follows. 


(3.53) 

In this case, the basis that diagonalize fermion mass matrices do not accord with those of 
the Yukawa couplings of the Higgs fields, thus flavor-changing effects are induced at radiative 
levels. The Higgs LEV effects contribute to the jU —> ey decay through the diagrams called 
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Barr-Zee type. They are relevant in case where masses of superpartilcles are heavy, while the 
CP-odd Higgs mass is relatively light. In the main part of this dissertation, since we assume 
that only the Bino and the sleptons are light, while the other superparticles are decoupled, the 
Higgs-mediated FCNCs are suppressed. 


3.4.2 Formulae in electron EDM 


The MSSM generically contains complex parameters, which cause CPV Especially, CP phases 
in lepton sector contribute to the electric dipole moment (EDM) of the leptons. The EDM can 
be calculated from the effective Lagrangian (2.6). The coefficient ^ 3 ( 0 ) can be evaluated in 
the same way as the muon g — 2. The one loop formulae for the electron EDM are given by 


^ ^ jji~± 


4871 ^ mi 

X=1 A=1 Vx 


1 TTlyO - 

EE 


L(e)* 

lAX 


9671" 


x=lA=1 


L(e)* 

lAX 



(3.54) 


(3.55) 


where Eqs. (3.54) and (3.55) come from the chargino-sneutrino and the neutralino-slepton 
contributions, respectively. 

Then, we discuss the CP violating phases without Lpy i.e., the slepton squared masses 
m~ are universal among the lepton flavor. In this case, the masses m| must be real because 
of hermitically of Lagrangian. Physical observables depend on CP violating phases of the 
following combinations [92,93]^ 

(t)i = arg = arg , (3.56) 

where M; (i = 1,2,3) are gaugino masses, A^ (/ = f, q) and BjU are coefficients of A- and B- 
terms, respectively. Since these new phases do not generically vanish, and then lead to large 
EDMs, they are strongly limited by experiments. Although it is difficult to suppressed these 
phases in general SUSY models, we assume that they are vanished by some mechanisms, for 
simplicity. 

Even if ~ 0 are realized, off-diagonal components of the lepton Yukawa matrices, 

which are transformed those of slepton mass matrices in mass eigenstates of leptons, could be 
new sources of CPVs, They are not suppressed if the slepton mass matrices are non-universal 
among the flavors because the super GIM mechanism does not work. In Sec. 4.2.3, we will 
discuss sensitivities of these phases to by investigating constraints of the electron EDM. 


3.5 Higgs Coupling in the MSSM 

If the muon g — 2 anomaly is due to the SUSY contributions, masses of smuons, electroweak 
gauginos, and higgsinos are required to be order of 100 GeV. Eurther, when the slepton mass 

^ Strictly speaking, there are other conbinations, e.g., arg(M;M*) (i y )). 
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matrices are nearly universal among the flavors, staus should be light. Since staus have rel¬ 
atively large Yukawa coupling to the Higgs boson, the stau contribution to Higgs coupling 
to di-photon becomes sizable when the left-right mixing parameter of the staus is large. In 
this section, we provide formulae for SUSY contributions to the Higgs coupling to di-photon. 
Detailed analysis will be carried out in Chapter. 4. 

The Higgs coupling to di-photon is composed of the contributions from SM particles and 
the MSSM. Theoretically, fc^ is represented as 


_ |^^^(SM) + ^^^(MSSM)| 
where is related to the Higgs decay rate as 


m 


rr) = 


2 2 
a 

102471 ^ 



(3.57) 


(3.58) 


The right-hand side in Eq. (3.57) is dominated by one-loop contributions. The SM contri¬ 
butions .^^^(SM) are composed of the electroweak gauge boson and the fermions at leading 
order. They are evaluated as [94] 


^^^(SM) = 


ShWW 

5. 


m 





c./Q?4 



(3.59) 


where the first term of the right hand side represents the contributions of the electroweak 
gauge boson, and the second terms are from fermions, m^, is the lightest (or SM-like) Higgs 
mass, N^ j and are the number of colors and the electric charge of particle /, respectively. 

The loop functions and aI\ are given by 

2 


A\{x) = 2-I-3x + 3x(2 — x)/(x), 

A^ix) = -2x(l-h(l-x))/(x), 

2 

where the function /(x) is 
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(3.61) 

(3.62) 


In X —> oo limit, the loop functions approach to 


—3/4, respectively. The 
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couplings ghww are written as® 


_ 

ShWW — 


“ 7^’ 


(3.64) 


Numerically, the W boson contribution is leading one in the SM. Among SM fermions, the top 
quark contribution is dominant because of large top Yukawa coupling. 

The MSSM contributions .^^^(MSSM) are composed of the sfermion ^^^(/), charged 
Higgs and chargino contributions respectively at the one-loop level. They 

are evaluated as [94] 





2g„ 


'^A ^A 


AA 



(3.65) 


The parameters mj, m^±, are masses of the sfermions, charged Higgs bosons, and charginos, 
respectively. The loop function is given by 

Aq(x) = x(l — x/(x)). (3.66) 

It approaches to A^ —> —1/3 in x —> oo limit. The couplings g^jj, ghn+H-’ Shx^x- 

g,f 7 = —{5m~ + 5m~ )±—{5m~ — 5m~ )cos20f ±5m~ sin267, 

6h/i.2/i.2 2 hi /rr 2 III /rr ^ fLR D 

f mzCos2p A 

gm+H- = ^ sin(^ - a) + —-— sin(a + /3) , 

2 cos ) 

Shx^xL = [^(-(OL)ii(OR)2iSina+(Oi);2(0^)i;Cosa), (3.67) 

where are the unitary matrices, which diagonalize the chargino mass matrix, as seen in 
Sec. 3.3.1. 5m~ and 5m~ are defined as 

/lL,RR flR 

Sml =-(m^+D7^J, 5ml =-ml , (3.68) 

® In the MSSM, there are additional mixing angle from the Higgs sector. The Higgs couplings to the SM 

particles in the MSSM are expressed by 

gv ruf —sin a 

g„„ = ^sl„(p-a). a„=__, (3.63) 

where a is mixing angle between the real part of the Higgs fields. In this dissertation, we assume that the CP-odd 
Higgs A is decoupled. In this case, the Higgs couplings to SM particles approach to the SM values [3 tt -I- 
Therefore, we can neglect the MSSM contributions to the Higgs couplings to SM particles. 
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r 

Figure 3.6: Feynman diagram of the stau contribution to the Higgs coupling to di-photon. 


where are D-terms, Dj = cos 2^ sin^ 6^j, and m~ is the left-right mixing 

of the sfermion mass matrices as 


m% 


flR 


—(m| — m|)sin20j. 


fi 


IT 


Here, 6j is mixing angles of the sfermions, defined as 


U7 


COS 6j sin 6j 
— sin 9j cos 


(3.69) 


(3.70) 


The mass matrix of the sfermion is diagonalize by Uj, as = diag(m~,m~). Here, 

^ ^ f f fi fi 

we neglect effects of flavor/CP violations. These contributions are expected to be suppressed 
since they are strongly constrained by lepton flavor experiments. 

The contributions of the charged Higgs bosons and the charginos becomes relevant when 
these masses are order of 100 GeV. For charged Higgs, the contributions are negligible since 
we assume that the CP-odd Higgs mass is decoupled. This assumption is reasonable since 
none of them is discovered. On the other hand, mass regions at several hundred GeV for 
charginos are still remained. As we will discuss in Sec. 4.3, the chargino contributions to fc^ 
might become In this dissertation, we take them as theoretical uncertainty. Detailed 

estimation will be discussed in Sec. 4.3. 

Then, let us consider the contributions of sfermions. Light sfermions with large left-right 
mixing, i.e., the stop, sbottom, and stau, can provide large contributions to the Higgs coupling 
to di-photon [95]. We assume that the colored superparticles are heavy enough to evade cur¬ 
rent LHC bound. In our assumption, stau may have a large contribution to the Higgs coupling. 
Fig. 3.6 shows the leading stau contribution to Higgs coupling to di-photon amplitude in large 
left-right mixing region. In this case, the leading stau corection is estimated as [96] 



(3.71) 



56 


CHAPTER 3. SUPERSYMMETRIC STANDARD MODEL 


where this approximation is applicable when the left-right mixing of the stau are extremely 
large (or m~^ » off-diagonal component of the stau mass 

matrices is comparable to the diagonal parts, ^(10)% correction to 5k^ can be realized. 

However, extremely large left-right mixing might suffer from vacuum instability [97,98]. 
As the left-right mixing increases, new charge-breaking minima become deeper, and our elec- 
troweak vacuum could decay them. By requiring that the lifetime of our vacuum should be 
longer than the age of the Universe, the left-right mixing is bounded from above [99,100]. 
Hence, the stau contribution to the Higgs coupling to di-photon is also limited. The upper 
limit is studied by several theoretical groups [96, 100, 101], and becomes about 15%. 

In Chapter 4, we will study allowed stau mass region by the Higgs coupling to di-photon 
K^, which is limited by the vacuum stability condition of stau-Higgs potential. Then, future 
prospects for the stau based on above mass region are discussed. 


3.6 Current Status 


Supersymmetry which provide a solution for the hierarchy problem, is required that the 
masses of superparticles are not so far above the electroweak scale, and has long been the 
leading candidate for new physics. This has motivated searches of superparticles at colliders, 
such as the Large Hadron Collider (LHC). Particularly, the main searches has been for strong- 
interactiong superparticles (colored superparticles), such as gluinos and squarks, because of 
large production cross section. They are produced through pp —> gg, gq, qq. Although these 
limit depend on the decays, masses (m^, is strongly limited by current LHC data, in case 
assuming some of simplified models (or decay modes). 

One of the benchmark models is the framework of minimal supergravity (mSUGRA), also 
knew as the constrained MSSM (CMSSM). It has 5 parameters, a unified scalar mass ttIq, 
a unified gaugino mass ^ unified trilinear scalar coupling Ag, the ratio of Higgs VEVs 

tan/3, and one discrete choice, the sign of p. Fig. 3.7 shows summary on constraints in 
mSUGRA/CMSSM. This results imply > 1.4TeV and > 1.6TeV [102]. 

Current mass bound of gluinos and squarks, ^ > 1 TeV, is consistent with the Higgs 
mass of 126 GeV. In the MSSM, the Higgs mass at tree level is about the Z boson mass, 
< rri^ cos^2p, because the quartic terms in the Higgs potential is dominated by gauge 
couplings. This contradicts the LEP bound [103], > 114.4 GeV. 

The relation to the Higgs mass receives quantum corrections. Relevant contributions come 
from top and stop loops. If the corrections are comparable to the tree level contribution, the 
Higgs mass of about 126 GeV can be realized. At one-loop order and taking in the decoupling 
limit (m^ » m^), the Higgs mass is evaluated as [10-14] 
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Figure 3.7: Constraints on mSUGRA/CMSSM from ATLAS collaboration at LHC with if' 
20fb“^ and -v/s = 8 TeV. Aq = —2mo, tan^ = 30, /r > 0 are taken to realize the Higgs mass 
L26GeV. In this figure, (mg, m^/ 2 ) plane of mSUGRA is limited from the jets + $j-. This 
figure is from Ref. [L02]. 
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where Mg = and Xf = Aj — jU cot /3. 

From Eq. (3.72),® it is found that there are two possibilities to enhance the Higgs boson 
mass as follows, 

• Ms » L TeV, Xf 0 (Heavy scalar top) [L0-L3] 

• Ms ^ LTeV, Xf VSMg (Maximal mixing) [L4] 


^ In fact, in order to evaluate the Higgs mass precisely, it is also important to take into account higher order 
corrections. However, qualitative consequence does not change. 
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Anyway, Mg > 1 TeV is favored to explain the Higgs mass, and is consistent with current mass 
bound of colored superparticles. 

These two phenomenological requirements, (A) the mass bound of colored superparticles 
is > ITeV, (B) the Higgs boson mass is 126 GeV, respectively, naively conflict with (C) the 
indication of the muon g — 2 anomaly where superparticles have a mass of ^(100) GeV, as 
seen in Sec. 3.3. In fact, some of the representative SUSY-breaking mediation mechanism such 
as mSUGRA and minimal GMSB cannot satisfy (A)-(G) simultaneously [15]. There are two 
possibilities to solve above inconsistency as follows, 

(i) Model building (extending MSSM). 

(ii) Model-independent approach. Particularly, we consider a split mass spectrum of the 
superparticles as 


ntg, ntq » mj, mg, m^, (3.73) 

where the squarks are heavy enough to explain the Higgs boson mass and satisfy the current 
LHG limit, while the non-colored superparticles are light to solve the muon g — 2 anomaly. 

In the case (i), some theoretical groups have proposed extended models to explain the 
muon g — 2 and the Higgs mass simultaneously. For example, an extension of MSSM with 
vector-like supermultiplets^° was considered to increase the Higgs boson mass by extra contri¬ 
butions [107-111], while explaining the discrepancy to the muon g — 2 [112,113]. Then, the 
collider phenomenology were studied in Ref. [114]. 

The scenario (ii) forces us to change strategy of superparticle search. The standard SUSY 
searches are based on productions of colored superparticles as mentioned in the beginning of 
this section. In this case, they are not always promising because the colored superparticles 
could be away from the LHG reach. On the contrary, the searches for weak interacting super¬ 
particles (non-colored superparticles), such as electroweak gauginos, Higgsinos, and sleptons, 
become important. If the muon g —2 anomaly is true, non-colored superparticles couple to the 
muon should have masses of ^(100) GeV, and are expected to be produced at LHG, through 
direct productions such as pp —> 

Now, let us review the searches for electroweak gauginos (chagrinos and neutralinos). 
The ATLAS [117,118] and the GMS [119] collaborations respectively reported results of their 
searches for above direct productions. The mass bound of the electroweak gauginos is shown 
in Fig. 3.8. The left panel of Fig. 3.8 is result of ATLAS collaboration, the right is from GMS. 
Although the results strongly depend on the mass spectrum of non-colored superparticles, the 
mass regions for the superparticles are beginning to be limited by current LHG data. 

LHG phenomenology of the superparticles which are relevant for the muon g — 2 has been 
also studied in Ref. [4]. Fig. 3.9 shows the LHG exclusion limits with the blue dotted lines 

Other extensions are additional U(l)' gauge S 5 nnmetry imposed on the Higgs fields [15], a large LR mixing 
of the stops by introducing a Higgs-messenger mixing [104,105], splitting F-terms between colored and non¬ 
colored messengers [106], and so on. 

This is previous our study. For detailed discussion, see Ref. [4]. 
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Figure 3.8: Mass bound of the electroweak gauginos, which are produced through pp —> 
~±,o~±,o^ by the ATLAS (left) and the CMS (right) collaborations, respectively. These figures 
are from Refs. [115,116]. 

at 95% CL, where the SUSY contributions to the muon g — 2 are dominated by chargino-muon 
sneutrino diagrams. When the slepton are heavier than the Bino, while are lighter than the 
Wino, a signal of 3 leptons plus large missing transverse energy is enhanced. Thus, it is found 
that the muon g — 2 regions are partly excluded by current LHC data, and it has been shown 
that almost all parameter regions are expected to be probed in near future by LHC at 14 TeV. 

However, the searches in Ref. [4] rely on the assumption that Wino is light. This assump¬ 
tion is not always necessary to explain the muon g — 2 discrepancy, when the SUSY contribu¬ 
tions is mainly from neutralino-smuon diagrams in Fig. 3.3 of Sec. 3.3. Among the neutralino 
contributions, Bino-smuon diagram as Fig. 3.4 (b) of Sec. 3.3 is proportional to the left-right 
mixing of the slepton. When only the Bino and the sleptons are light, and the left-right mixing 
is large, the Bino-smuon contribution is enhanced. Eventually, we notice that the models in 
which only the Bino and the sleptons are light, are "minimal" SUSY models to solve the muon 
g — 2 anomaly. 

If the Wino is heavy, collider searches differ significantly from those in Ref. [4], thus it 
is important to investigate the phenomenology of the "minimal" SUSY models to explain the 
muon g — 2 discrepancy. In the next chapter, we will study phenomenology of this "minimal" 
models, where only the Bino and the sleptons are light. 
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M2[GeV] 


Figure 3.9: Current LHC bounds on the SUSY g — 2 explanations [4], by the production 
PP X 1 X 2 + $T the data of f L/f = 13.0fb“^ [120]. The parameters are taken 

jU = M 2 = 2Mi, mj^ = 3TeV, and tan p = 40, respectively. In the orange (yellow) regions, the 
muon g — 2 discrepancy is explained at Icr (,2a) level, where Icr (2cr) mean that the SUSY 
contribution to the muon g — 2 are (26.1 i 8.0) x 10“^° ((26.1 i 16.0) x 10“^°), as will see 
in Sec. 4.1.2. The blue dotted line represents the LHC exclusion limit at 95% CL. Including 
30% errors of Parton Distribution Function (PDF) and scale uncertainties, the limit distribute 
in the blue hatched regions. The LSP is Xi in the regions above the black thick lines, while 
the sneutrino is lightest below them. Sleptons are lighter than Xi (X 2 ^ below the black dotted 
(dashed) lines, respectively. We find that the muon g —2 regions are partly excluded by current 
LHC data. 








Chapter 4 

Prospects for Slepton Searches 


This chapter is the main part of this dissertation and is based on [1-3]. In Sec. 4.1, we 
introduce an effective theory in which only the Bino and the sleptons are light, while the 
other superparticles are decoupled. In this setup, the left-right mixing becomes large in order 
to solve the muon g — 2 anomaly. Such a large mixing spoils the stability the electroweak 
symmetry breaking vacuum. We also discuss the vacuum meta-stability condition of the Higg- 
slepton potential and the slepton mass bound. 

Searches for sleptons depend on the mass spectrum of the sleptons. We will show the left- 
right mixing of the slepton is constrained by the vacuum meta-stability condition of the slepton 
in Sec. 4.1.3. Then, it is found that the slepton mass is bounded from above in order to explain 
the muon g — 2. Since the most severe constraint is obtained from the stau stability condition 
because of large tau Yukawa coupling, we consider the two following mass spectrum, 

(1) Universal mass spectrum (Sec. 4.2.2, 4.3) 

^e = ^n = (4-1) 

(2) Non-universal mass spectrum (Sec. 4.2.3) 

m-^ = mj2 (4.2) 

In the Universal case (1), we will show that smuons are limited to be < 500GeV 
to solve the muon g — 2 anomaly. This is within kinematical reach of the LHC and/or ILC 
sensitivity. The searches is discussed in Sec. 4.2.2. Furthermore, light staus also affect the 
Higgs coupling to di-photon k^, whose sensitivity is expected to be improved by HL-LHC and 
ILC. Once the excess from the SM prediction of is observed, we can estimate not only the 
stau contribution to k^, but also their masses and mixing. In Sec. 4.3, we discuss the stau 
searches in ILC. 

On the other hand, the vacuum condition is relaxed when staus are much larger than other 
sleptons. In this case, smuon masses of ^(l)TeV are allowed to solve the muon g —2 anomaly, 
and outside the kinematical reach of HL-LHC and ILC. However, such a non-universal mass 
spectrum generically induces sizable lepton flavor violation and the CP violation. In Sec. 4.2.3, 
we discuss the sensitivity to the future lepton flavor experiments. 
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4.1 Foundation 

4.1.1 Setup 

In this chapter, we consider a low-energy effective theory, in which only the following super¬ 
particles are light. 


B, h, 4- (4.3) 

Here, i denotes the selectron and the smuon (and the stau, depending on the mass spectrum). 
The model parameters are as follows. 

Ml, m? , m? , m? . (4.4) 

’^R ’^LR 

Here, Mj is the Bino mass. On the other hand, m? and m? are soft SUSY-breaking masses 

^R 

of the left- and right-handed sleptons, respectively, m? is off-diagonal components of the 

^LR 

slepton mass matrices. This is the "minimal" SUSY model to solve the muon g — 2 anomaly. 

All colored superparticles are set to be very heavy. In fact, none of them have been 
discovered at LHC. The Higgs boson mass of 126 GeV favors the scalar top masses to be 
^(10 — 100)TeV, if the trilinear coupling of the top squark is suppressed. Similarly, the heavy 
Higgs bosons of the two Higgs doublets are assumed to be heavy for simplicity. In this disser¬ 
tation, all of them are considered to be decoupled, although the results of this dissertation are 
almost independent on their mass spectrum.^ 

In this dissertation, we consider the case that the Higgsino mass is large. The SUSY contri¬ 
butions to the muon g — 2 are dominated by the Wino-Higgsino-muon sneutrino and/or the 
Bino-smuon contributions. In our setup, the Bino-smuon contribution is dominant in order 
to explain the discrepancy of the muon g — 2. It is enhanced by the left-right mixing of the 
smuon, which is determined by the muon Yukawa coupling, the Higgsino mass jji, and tan/3. 
Since too large tan ^ spoils perturbativity of the down-type Yukawa interaction, jji is favored 
to be large. Therefore, we focus on the large Higgsino mass region. 

When M 2 and jjl > mj^, the chargino-muon sneutrino contribution to the muon g — 2 is 
proportional to l/M 2 /r, where M 2 is the Wino mass. Thus, it is suppressed by large (M 2 , ju) 
case. The Winos are also supposed to be decoupled for simplicity.^ 

Before proceeding, let us comment on the left-right mixing of the slepton. It includes the 
scalar trilinear coupling of the slepton, A^, as well as jU and tan^. The Bino-smuon con¬ 
tribution to the muon g — 2 could be enhanced by A^ with tan /3 kept small. However, 

^ The non-decoupling effect of the heavy superparticles appears the Bino coupling with the smuons as will 
see in Sec. 4.1.2. They yield correction of 5-10% to the SUSY one-loop contribution to the muon g — 2. In this 
dissertation, masses of the heavy superparticles are set to 30 TeV. 

^ Note that the setup of this dissertation is meaningful when the Wino is decoupled. When M 2 < 1 TeV, 
the chargino-muon sneutrino contribution can be effective, and cannot be neglected. Therefore, the analysis 
dramatically changes due to large theoretical uncertainty in this case [4] . In this dissertation, we set to M 2 = 
30 TeV in order to suppress the chargino-muon sneutrino contribution. 
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Figure 4.1: The Bino-smuon contribution to the magnetic dipole operators. 


this requires to be extraordinary large. If the trilinear coupling is universal among the 
matter scalar fermions, this implies extremely large trilinear coupling for the stop sector, re¬ 
sulting in either too large Higgs boson mass or rapid decays of our electroweak vacuum into 
charge/color breaking vacua. In this dissertation, A^ is set to be zero for simplicity, and the 
left-right mixing of the smuon is determined by jji and tan^. 

4.1.2 Muon g — 2 

In this section, we discuss the SUSY contribution to the muon g — 2 in the minimal SUSY 
model. The feynman diagram for the Bino-smuon contribution to the muon g — 2 is shown in 
Fig. 4.1, where the lepton £ = ju. In the mass insertion approximation, it is estimated as 





(lOOGeV)^ 


ml ml /M.ix 

Mr' 


at the leading order of tan The loop function f^ix, y) is defined as Eqs. (3.44). It is noticed 
that the contribution is proportional to the left-right mixing of the slepton, ml . 


LR 



(4.6) 


where is the vacuum expectation value (VEV) of the up-type Higgs field and tan ^ is 
the ratio of Higgs VEVs, tan^ = {Hj/{H^). The parameter is a correction to the lepton 
Yukawa coupling constant [121]. It appears when the lepton Yukawa coupling in the MSSM 
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is matched to the physical lepton mass or the lepton Yukawa coupling in the SM. In the 
low-energy effective theory, it is evaluated as 

Af tan ^ • I m|J , (4.7) 


where the diagrams including Higgsino propagator are neglected, because they are suppressed 
in large ju regions. Also, no-tan^ enhanced terms are discarded. The loop function I {a, b,c) 
is defined as 


/(a, b,c) = - 


ab log (f) + he log + ca log 
(a — b)(b — c)(c — a) 


(4.8) 


In particular, when jUtan/3 is very large, becomes as large as or larger than ^(0.1 — 1). 

The correction denotes leading contributions of higher order corrections. It is esti¬ 
mated as 


1 + = 



4a 

n 



2ay Ab H—a, 
4 ^ 


In- 


M, 


soft 


^soft 


(4.9) 


In the right-hand side, the first bracket is QED corrections to the muon g — 2 [122]. They 
can be estimated by renormalization group contributions the magnetic dipole operator from 
the (slepton or Bino) soft mass scale rrig^f^ to the muon mass scale. In the numerical analysis, 
is chosen to be the smuon mass, and this correction is ~ 10%. Non-logarithmic terms 
evaluated in Ref. [122] are found to be very small in our parameter regions. 

The second bracket in Eq.(4.9) is corrections to the Bino couplings with the smuons.^ 
When SUSY is unbroken, the gaugino coupling is equal to the gauge coupling constant thanks 
to SUSY. This equality is violated after heavy superparticles are decoupled at a scale The 
Bino-muon-smuon interactions are 

-^int Bpi K + + h.c., (4.10) 

where the coefficients are the Bino couplings as 


gL=gY + SgL ^ gr 
gR = gY + SgR ^ gy 


1 H-I ayAb H— 

4n 


1 I A U 1 ^Soft 

1 H-Ab In- 

471 


In 


■^soft 
^soft _ 


(4.11) 

(4.12) 


Eor derivation, see Appendix B. Here, the terms of ay are corrections to the Bino self-energy 
which are SUSY analog of the oblique corrections [124-128]. They are called the super¬ 
oblique corrections. The U(l)y gauge coupling constant, gy, is evaluated at The cor¬ 

rections are represented by the difference between the beta functions of the U(l)y gauge and 

^ Ref. [123] also studied non-decoupling two-loop contributions from sfermions, which partially overlaps with 
our discussion. Our results are consistent with those of Ref [123]. 
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Bino couplings. In the minimal model, SM particles and light sleptons contribute to the beta 
functions. Thus, the coefficient becomes Ab = 41/6 — where is number of 

the generations of light sleptons; for instance, n^iepton = 3 if all the sleptons are light, and 
rigigpion = 2 when the staus are decoupled. On the other hand, the term of a 2 is non-oblique 
corrections after the Wino decoupled. In case that the gluinos are decoupled, similar calcu¬ 
lation is evaluated in Ref. [129]. The weak boson couples only to the left-handed (s)leptons 
at the Bino-muon-smuon vertices. Since the Bino-smuon contribution to the muon g — 2 in¬ 
cludes both and g^^, the loop correction (4.9) is obtained. It yields a correction of 5 — 10% 
for Msoff = 10 — 100 TeV with ~ 100 GeV. In our analysis, only the logarithmic terms are 

considered. Non-logarithmic terms are expected to be suppressed. 

In Eq. (4.9), we include the corrections that can be as large as or larger than ^(10)%. 
Other 2-loop corrections are unknown or expected to be small. The SUSY corrections to the 
SM one-loop diagram is calculated in Refs. [130-134]. When the superparticles except those 
of (4.4) are decoupled, these contributions are negligibly small. Other corrections that have 
not yet been estimated include electroweak and SUSY two-loop contributions to SUSY one- 
loop diagrams. They might provide 5a^ ~ 10“^°, according to Ref. [135]. In addition, non- 
logarithmic corrections to could be a few percents of the SUSY one-loop contributions, 
similarly to the discussions in Ref. [135]. 

We comment on other one-loop contributions to the muon g — 2. The Bino-Higgsino- 
smuon contributions (3.41), (3.42) can be < ^(10“^°) when the Higgsinos are relatively light 
due to the vacuum stability bound of the stau-Higgs potential in Sec. 4.1.3. We include these 
contributions in the numerical analysis for completeness.'^ On the other hand, the chargino- 
muon sneutrino contributions (3.40), which are dominated in the case of nearly degenerate 
mass spectrum, are less than ^(10“^^) for M 2 > lOTeV, i.e., negligible. We do not include 
them in the numerical analysis. 

Summarizing the above, the theoretical uncertainty of the SUSY contribution to the muon 
g — 2 can be estimated at ^(10)%. It is smaller than the current experimental and theoretical 
(SM) uncertainties. Here and hereafter, we understand that the results of this dissertation 
have uncertainty of ^(10)%, which depends on the heavy mass scale and unknown 

higher-order corrections. In this dissertation, we set to = 30 TeV, for simplicity. 

In Fig. 4.2, we show contours of the SUSY contributions to the muon g — 2. The horizontal 
and vertical axises are the lightest smuon mass and jU, respectively. The parameters are 
set as tan^ = 40 and = 30 TeV. In the right panel of Fig. 4.2, a part 

of the parameter region of the left panel is showed. In the orange (yellow) regions, the SUSY 
contribution to the muon g — 2 becomes (26.1 i 8.0) x 10“^° ((26.1 i 16.0) x 10“^°). In this 
dissertation, we define these regions as Icr (2cr) regions. Hereafter, we will use the above 
definition without notice. From Fig. 4.2, it is found that they are enhanced by large jj,, and the 
smuon masses can be 1 TeV for jU = ^(10 —100) TeV. This is contrasted to the chargino-muon 
sneutrino contributions to the muon g — 2, where jji is favored to be small [4] . On the other 

This contribution can dominate the SUSY contributions to the muon g —2, when is small while decoupling 
the Wino. Since they are enhanced only by tan^, superparticles are required to be light to explain (2.34). They 
are detectable in colliders. In particular, the Higgsino production can be significant. 
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Figure 4.2: The SUSY contributions to the muon g — 2 and the vacuum stability bounds are 
shown. In the orange (yellow) regions, the muon g — 2 discrepancy (2.34) is explained at 
the lc7 (2cr) level, where Icr (2cr) mean that the SUSY contribution to the muon g — 2 are 
(26.1 it 8.0) X 10“^° ((26.1 it 16.0) x 10“^°). The black solid lines represent the upper bound 
on jJL from the vacuum stability bound of the stau-Higgs potential. The stau masses are set to 
be m^/rriji = 1,2,5,10 and 20 from bottom to top. Above the black dashed line, becomes 
larger than unity. The parameters are tan ^ = 40 and = 30TeV. In the 

right panel, a part of the parameter region of the left panel is magnified. 


hand, detailed dependences on the superparticle mass spectrum are determined by the loop 
function (3.44) and the vacuum stability condition. They will be discussed in Sec. 4.1.3. 


4.1.3 Vacuum Stability 

As shown in Sec. 4.1.2, the Bino-smuon contribution to the muon g — 2 is enhanced by a large 
left-right mixing of the smuon. However, too large left-right mixing spoils the stability of the 
electroweak symmetry breaking vacuum (EW vacuum). In this section, we discuss the vacuum 
stability condition of the sleptons and slepton mass region using the condition. 

Let us consider the scalar potential for the neutral component of the up-type Higgs, H^, 
the left-handed slepton, and the right-handed slepton, It is given as [1] 


V = (m^^+/t)|Hj +m~J£J +m~J4| - + h.c.) + Yj£*4| 

+|-(|?J2 + \Hjy + ^(|Zj2 - 2|4P - \Hff + ^^^5 h\hX, 


8 


(4.13) 
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where is the soft SUSY-breaking mass of the up-type Higgs. When tan 6 and the CP-odd 
Higgs mass, rrij^ are large, the terms including the down-type Higgs, can be neglected.^ 
Further, only the renormalizable terms are taken into account in this dissertation. Higher¬ 
dimensional terms depend on the superparticles that decouple at which is mentioned in 
Sec. 4.1.2. Since is very large, their contributions to the vacuum meta-stability condition 
are considered to be small. 

We take account of the RG-improved term in the scalar potential, which is the last term in 
Eq. (4.13). The leading correction for 5^ is represented as 




2n^ g2 + ga 


log 




2 2 

t^ to 


mf 


(4.14) 


where is the top Yukawa coupling, and and are the top and stop masses, respectively. 
The Higgs mass is enhanced if is sizable. In our analysis, we set 5^ ~ 1 to reproduce the 
Higgs mass of 126 GeV. 

The coefficient of the quartic terms in the scalar potential, is somewhat important. 
At the tree level, = m^/(vcos^) (m^/v)tan/3 in a large tan^ limit. As mentioned in 
Eq. (4.6), It is modified at the one loop level as follows: 


rrif 


m( tan [3 


V cos 13 1 + V 1 -I- 

The dependence on Y^ to the vacuum stability condition will be discussed later. 
Expanding around the EW vacuum, H^ = v^ + -^h^, we have 


(4.15) 


3- Sv ~ ~ f tn~ 


(■LL 

2 


— m: 



32 


) 2 2 2 


{sV2vX+K) + 




8 


(4.16) 


where m~ are the squired masses of the left- and right-handed sleptons as seen in (3.20) of 

iLL,RR 

Sec. 3.3.1, and m~ is the left-right mixing of the sleptons defined in Eq. (4.6). 

Ur 

The scalar potential has the ordinary electroweak symmetry breaking minimum (EW vac¬ 
uum), = (-L — — 0- the trilinear coupling of the sleptons and the Higgs boson, m~ 

Ur 

^ Authors in Ref. [100] analyzed the scalar potential including the down-type Higgs field for completeness. As 
tan^ approaches to 1, the effects of the down-type Higgs become relevant. However, we will not consider such 
small tan [3 region since it is difficult to explain the muon g — 2 anomaly. 
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increases, disastrous charge-breaking minima (CB vacuum), 7 ^ 0 or 0 in the scalar po¬ 
tential become deeper, and our EW vacuum could decay them by quantum tunneling effects. 
By requiring that the lifetime of the EW vacuum should be longer than the age of the Universe, 
the mixing m~ is bounded from above. 

t-LR 

The vacuum transition rate from the false vacuum (EW vacuum) to the true vacuum (CB 
vacuum) can be evaluated by semiclassic technique [136,137]. In this technique, the imagi¬ 
nary part of energy of the false vacuum corresponds to the transition rate (see Appendix A). 
The vacuum transition rate per unit space-time volume is evaluated as follows, 

r/V=A-e-®, (4.17) 


where the prefactor A is the fourth power of the typical energy scale in the potential, and 
expected to be roughly (100 GeV)"^. Although it is hard to precisely calculate A, the decay rate 
is not sensitive to it. On the other hand, the Euclidean action which is evaluated on the bounce 
solution, B, is important to estimate the decay rate. It is defined as [99] 


B = S£[0(r)]-SE[<^^], (4.18) 

where S^ is Euclidean action, ^ is a bounce configuration, and (f)^ is field values at the false 
vacuum. When T/V is smaller than the fourth power of the present Hubble parameter Hq = 
1.5 X 10“"^^ GeV, (i.e. the lifetime of the false vacuum is longer than the age of the Universe), 
the false vacuum becomes unstable. The vacuum meta-stability condition is approximately 
given as B > 400 if we assume A ~ (100 GeV)"^. 

The public package CosmoTransitions 1.0.2 [138],^ which is the numerical package 
to analyze cosmological phase transitions for single and multiple scalar fields, is used for the 
numerical evaluation of the Euclidean action B. The scalar potential is analyzed at the zero 
temperature. The vacuum mata-stability conditions correspond to constraint of m~ . The 

^LR 

fitting formula of the condition is obtained as [ 1 , 2 ] 


m~ 


< 


Vi 


1.01 X 10^ GeV^mj^ mj^ + 1.01 x 10^GeV(m|^ + 1.03m^^) 


-2.27xlO^GeV^ + 


2.97 X lO^GeV^ 




1.14 X 10®GeV^ 


+ 


0.983 


m~ 


m~ 


(4.19) 


The fitting formula reproduces results of CosmoTransitions at better than the 1% level. 
Note that an asymmetry of the coefficients of and mj^ is due to the D-term potential, 
i.e. the terms including and g^ in Eq. (4.13), although the effects of this asymmetry are 
negligible in numerical analysis. 

Let us compare Eq. (4.19) to the result in Ref. [99]. The meta-stability condition of the 
stau-Higgs potential was evaluated as [99] 

|/xtan^| < 213.5- 17.0(m;f^ 


^ In order to evaluate the bounce solution and the B numerically, we used modules tunnelinglD.py and 
pathDef ormation.py, and set the parameter a = 3 in the modules. 
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Figure 4.3: tan/3 dependence of the vacuum stability condition of the stau-Higgs potential. 
Here, we set to m= = = 300 GeV. 



+ 45.2 X lO-^GeVHm^^ - - 1.30 x 10^ GeV, (4.20) 

by numerical calculation. It is valid in the regions of < 600GeV, and fails as either 

the left- or right-handed slepton mass becomes extremely large. The fitting formula (4.19) 
is not only consistent with Eq. (4.20) in light slepton mass regions, but also applicable when 
masses of the sleptons are ^(l)TeV. 

The fitting formula (4.19) is universal for all the slepton flavors up to corrections, 17 ^ ~ 1. 
Given the soft scalar masses, the meta-stability condition is independent of the flavor except 
through Yukawa interactions in the quartic terms of the scalar potential. Since the quartic 
terms are dominated by gauge interactions, 17 ^ changes little around unity. It also depends 
on tan/3 when the Yukawa coupling is large. In practice, 17 .^ = 1 is set for the stau with 
tan^/(l + A.J.) = 70 [1, 99]. Numerical result of 17 for the stau with various tan/3 is found 
in Fig. 4.3. Here, although the stau mass is taken to be = 300GeV, the result is 

(most) independent of the stau mass. For instance, it is rj.,. 0.94 for tan/3/(l-l-A.j.) = 40. On 

the other hand, we obtain 0.88 for the smuons and the selectrons, which is independent 
of the flavor and tan ^ because of small Yukawa couplings [2] . 

The theoretical uncertainty of the fitting formula (4.19) is expected to be at most ^(10)%. 
First, the meta-stability condition corresponding to B has only logarithmic dependence of Hq 
and A. Thus, even if the value of Hq or A change ^(10)%, the change of the condition is less 
than one percent. Then, thermal corrections can change the transition rate. They could be 
significant (but < ^(10)%) when the staus are light, for instance, for < 200 GeV [139, 
140] in the scalar potential for the SM Higgs and the lightest stau. Although the corrections 
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Figure 4.4: The SUSY contributions to the muon g — 2 as a function of the lightest smuon 
mass, and the lightest neutralino mass, m^o. In the orange (yellow) regions, the muon 
g — 2 discrepancy (2.34) is explained at the let (2cr) level. The left-right mixing is maximized 
under the vacuum stability condition. The parameters are set to mj^ = mj^, tan /3 = 40 and 
Mgoft = lOTeV. The stau soft masses are equal to those of the selectrons and smuons. The 
region below the green line is excluded by LHC. 


in case of three-dimensional potential are not discussed, it is expected that the results are 
almost unchanged, since the trilinear coupling of the lightest stau with the Higgs boson is 
more important for the vacuum meta-stability condition. For more precise evaluation, detailed 
analysis is needed. 

The most severe constraint on jUtan/3 is obtained from the stability condition of the stau- 
Higgs potential. This is because the left-right mixing is proportional to the Yukawa coupling. 
In Fig. 4.2, we show the upper limits of jutan ^ for tan (3 = 40. Combined with the muon g — 2, 
the smuon masses are bounded from above for given stau masses. For = m^, the 

lightest smuon is limited to be < 300 (420) GeV at the Icr (2cr) level of the muon g — 2. 

The SUSY contributions to the muon g — 2 depend on superparticle mass spectra mainly 
through the loop function (3.44). In Fig. 4.4, we show contours of the muon g —2 as a function 
of smuon and neutra li no masses. In the orange (yellow) regions, the muon g —2 discrepancy is 
explained at the la (2a) level. The left-right mixing is maximized with satisfying the vacuum 
stability condition. Here, the stau soft mass is supposed to be degenerate with that of the 
smuon, nig = rriji = m^. It is found that the lightest smuon mass can be 330 (460) GeV at the 
la (2a) level of the muon g — 2, when the lightest neutralino mass, m^o, is almost a half of 

A 1 

Note that m^o is almost equal to M^. On the other hand, m^o ~ 560GeV is realized for 
250 GeV at the 2a level of the muon g — 2. 
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Figure 4.5: The left-right mixing of the smuon (blue solid) and the stau (red solid) as a 
function of jU. The blue (red) dashed line represents the vacuum meta-stability bound with 
p = 0.88 (0.94). Above the black dashed line, sleptons become tachyonic. Below the blue 
dotted line, < 1 is satisfied. The parameters are mj = mj = 2M^ = 500 GeV, tan^ = 40 
and Msoft = SOTeV. 


For a fixed value of the lightest smuon mass, the SUSY contributions to the muon g — 2 is 
maximized when the left-handed smuon mass, is equal to the right-handed smuon mass, 
Eq. (4.5) depends on the left- and right-handed smuon masses as fNi^,y)/xy with 
X = and y = rrAjM^. When xy is fixed, it is maximized for x = y and rapidly 

decreases for x ^ y. On the other hand, the vacuum stability condition (4.19) is relaxed for 
mj^ 7 ^ mj^. Since the former dependence is stronger than the latter, the SUSY contributions 
to the muon g — 2 become largest when is equal to for given and tan/3. In other 
words, the orange/yellow regions in Fig. 4.4, i.e., the parameter regions favored by the muon 
g — 2, shrink toward smaller superparticle masses for ^ 

The vacuum stability condition is relaxed as staus become heavy. When they are decoupled, 
the constraint almost disappear, and extremely large p, is allowed. As p increases, the left-right 
mixing (4.6) approaches to 


m~ =-ptanp 

°LR 1 -|- A^ 


471 1 

ay Ml I(Mt, m~ ,m~) 


(for p 


oo). 


(4.21) 


In Fig. 4.5, the left-right mixings of the smuon and the stau, \m~ |, are plotted as a function 

^LR 

of p, by the blue and red solid lines, respectively. Since m~ is proportional to the Yukawa 

Ur 

coupling, is larger than When the staus are light, p is bounded from above by 
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Figure 4.6: Mass bound for smuon when the staus are decoupled. In orange (yellow) regions, 
the muon g — 2 discrepancy is explained at Icr (2cr) level. The parameters are 
tan /3 = 40, and M^oft = 30 TeV. The left-right mixing is limited by the vacuum stability condi¬ 
tion of the smuon-Higgs potential. On the black line, = 1 is satisfied. 


the stau stability condition (red dashed line). For heavier staus, the red dashed line (as well 
as the stau tachyonic bound or the black dashed line) is lifted up, and larger p, is allowed. 

Even if the staus are decoupled, the left-right mixing is limited by the vacuum stability 
condition of the smuon-Higgs potential. Fig. 4.6 shows contours of the muon g — 2 under 
the smuon stability condition, . In the outside of the black-dashed line, the constraint (4.19) 
is more severe than Eq. (4.21), i.e., is satisfied for given smuon masses. 

We show that the lightest smuon mass can be as large as 1.4 (1.9) TeV at Icr (2cr) level of 
the muon g — 2. Since it is very weak, this bound does not appear in the parameter range of 
Eig. 4.2. 

In p ^ oo limit, the situation is essentially the same as the soft Yukawa coupling mod¬ 
els [141,142], where the Yukawa couplings at tree level equal to zero, and the fermion masses 
are generated by only radiative corrections. In this dissertation, the vacuum stability condi¬ 
tion is analyzed seriously and is found to restrict the SUSY contributions to the muon g — 2, 
resulting in tighter bound on the smuon mass than Ref. [142]. 

Note that is larger than unity when p is very large. Then, the radiative correction 
exceeds the tree level contribution to the lepton mass. Eor » 1, the Yukawa coupling 
constant is suppressed by A^, as observed by Eq. (4.6) . This may be disfavored by a naturalness 
argument [143]. In Eig. 4.2, A^ becomes larger than unity above the black dashed line, and 
above the blue dotted line in Eig. 4.5. 

So far, we discussed that the slepton masses are bounded from above by the vacuum meta- 



4.1. FOUNDATION 


73 


stability condition and the muon g — 2. Furthermore, when the slepton masses are almost 
same among lepton flavor, the stau mass region might be also detected by the Higgs coupling 
to di-photon jc^. In the next subsection, we discuss the stau mass region. 


4.1.4 Stau Mass Region 

We defined the minimal SUSY model, where only the Bino and the sleptons are light, while 
the other superparticles are decoupled in Sec. 4.1.1. If sleptons have nearly degenerate mass 
spectrum, staus also become light and the stau contribution to the Higgs coupling to di-photon 
JCj, might be sizable. Maximal value of fc^ (or the decay rate of r(h —» yy)) by the vacuum meta¬ 
stability condition have been discussed by some theoretical groups [96, 100, 101]. In contrast, 
we discuss the stau mass region where the Higgs coupling is deviated from SM prediction 
in the viewpoint of the vacuum meta-stability then classify the stau contribution. Our study is 
expected to provide some hints whether the deviation is from the stau contribution or not, if 
an excess of fc^ is observed. 

In this section, we discuss the stau mass region based on the Higgs coupling and the 
vacuum meta-stability condition of the stau-Higgs potential. The staus are characterized by 
the three parameters as follows. 


m^f.m 


'^2 


e. 


(4.22) 


where the left-right mixing angle is associated with the off-diagonal element of the stau 
mass matrix as (cf. Sec. 3.5) 



(4.23) 


Here, is assumed. The stau contribution to the Higgs coupling is determined, once 

the parameters (4.22) are given. They are constrained by the vacuum meta-stability condition. 

In Fig. 4.7, we show contours of = ?c^ — 1 as a function of 6^ and where the 
green solid lines are contours of for given The stau contribution depends on 9^ and 
is maximized when sm26if is close to unity (9^ ~ 7t/4) for fixed and Also, 5fc^ is 
enhanced by larger On the other hand, if T 2 is extremely heavy (T 2 = ^(l)TeV), the stau 
contribution to becomes insensitive to and controlled by and 9^. 

In Fig. 4.7, the red regions are excluded by the vacuum meta-stability condition. Eq. (4.19) 
gives an upper bound on for given and Combined with Eq. (4.23), the mixing 
angle 9^ is constrained as a function of and When is small, the angle is not 
limited by the vacuum meta-stability condition, and 5k^ is maximized when sin 20^ = 1 is 
satisfied. It is found that 5k^ becomes largest just below the red region with sin 20^ = 1 in 
each panel of Eig. 4.7. On the other hand, the vacuum meta-stability condition constrains 
the stau mixing angle for large The maximal value of decreases, as increases. 
When T 2 is extremely heavy, the vacuum meta-stability condition becomes insensitive to 
and determined by . This is because, in the decoupling limit, T 2 does not contribute to the 
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Figure 4.7: Contours of 5k^ are shown by the green solid lines. The lightest stau mass is 
90 GeV in the left panel and 150 GeV in the right. Here, we set to tan 13 — 20 and A^ = 0. The 
red regions are excluded by the vacuum meta-stability condition (4.19). The blue regions are 
excluded by the chargino search at LEP experiment. 


field configuration of the bounce solution to derive the vacuum meta-stability condition. The 
maximal value of 5 k^. is determined by . 

In the analysis, we choose tan^ = 20, A^. = 0 and M 2 = 500 GeV, where M 2 is the Wino 
mass.^ The stau contribution to and the vacuum meta-stability condition are almost inde¬ 
pendent of them, once is given. Rather, they are included in the definition of in 
association with the Higgsino mass parameter jU. For fixed jJL becomes smaller, as tan (3 
increases. When the charginos are light, they can affect the Higgs coupling [95,144]. Their 
contribution to is taken into account for completeness. It is at most a few percents in the 
vicinity of the blue region and much less than 1% around the red region in Fig. 4.7. The blue 
region is already excluded by LEP experiment [145], where the lightest chargino mass is less 
than 104 GeV. 

Next, let us study the stau mass region. In Fig. 4.8, we show contours of 5k^ by the green 
solid lines. Here, 5k^ is maximized with satisfying the vacuum meta-stability condition (4.19) 
at each (m^^,m^^). Each contour is composed of the two regions. In the left region of the 
peak, i.e. small region, sin20;f = 1 is satisfied. The stau contribution to is enhanced 

^ We consider the chargino contribution as theoretical error in this analysis. Detailed estimation for the error 
will be discussed in Sec. 4.3.2. 
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Figure 4.8: Contours of 5k.^ are shown by the green solid lines. At each point, is 
maximized under the vacuum meta-stability condition. Here, tan^ = 20, = 0 and 

7i/ 4 < < 7r/2 are taken. The blue dashed lines are contours of sin20;f. In the left re¬ 
gion of the blue dotted line (the leftmost blue line), = 1 is satisfied. 


when is larger, as explained above. On the other hand, in the right region of the peak, 0^ 
is limited by the vacuum meta-stability condition. Here, sin20;f is less than unity. Contours 
of sin20;f are shown by the blue dashed lines in Fig. 4.8. As already found in Fig. 4.7, fc^ is 
enhanced, when is smaller. 

In the Figs. 4.8, we also choose tan^ = 20, A^. = 0 and M 2 = 500 GeV. The results are 
almost independent of them except for the region in the vicinity of . When is 

very close to the Higgsinos become light, because the stau left-right mixing parameter 
tends to be small (see Eq.(4.23)). Then, the charginos can contribute to Otherwise, their 
contribution is negligible in Fig. 4.8. In the figure, it is also supposed that the lightest stau is 
mainly composed of the right-handed component, n/4 < 6^ < 7r/2. As mentioned above, the 
stau mass region in Fig. 4.8 is almost insensitive to this choice. 

Currently, the measured values of at LHC are consistent with the SM prediction. The 
uncertainties are 15% (ATLAS) [72] and 25% (CMS) [73]. As found in Fig. 4.8, they are not 
precise enough to probe the stau contribution for > 100 GeV. In future, the sensitivity 
will be improved very well, as mentioned in Sec. 2.5. It is expected that LHC measures at 
about 7% and 5% for the luminosities, 300fb“^ and 3000fb“\ respectively with ^/s = 14TeV 
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[146,147]. If the measurement of Br(h —> 77 )/Br(h —> ZZ*) at HL-LHC is combined with the 
measurements of the Higgs couplings at ILC, it was argued that the uncertainty of can be 
reduced to be about 2% (Icr) at 250GeV ILC with ^ — 250fb“^ [74]. If the luminosity is 
accumulated up to 2500fb“^ at 1 TeV ILC, it has been estimated that the accuracy of can be 
better than 1% [74]. 

It is noteworthy that, once an excess of is measured, the mass region of staus are 
determined from Fig. 4.8. From the joint analysis of 250 GeV ILC and HL-LHC, Sk^. is expected 
to be measured with the uncertainty of 2% at the Icj level. If 5k^ is measured to be larger 
than 4%, the upper bound is obtained as < 200 GeV.® Such a stau can be discovered 
at 500 GeV ILC. In fact, the stau is detectable up to 230 GeV at ILC with = 500 GeV and 
^ = 500fb“^ [148]. On the other hand, if 5k^ is measured to be 2% (1%), the stau mass is 
predicted to be less than 290GeV (460GeV). This is within the kinematical reach of 1 TeV ILC. 
Therefore, if the stau contribution to is large enough to be measurable, the stau is predicted 
to be discovered at ILC.^ 

Let us mention the case when the heaviest stau is very heavy. In contrast to the lightest stau 
T^, the heaviest stau T 2 can be decoupled with enhanced and the vacuum meta-stability 
condition satisfied. In Fig. 4.8, 5k^ is insensitive to and determined by for very large 
In the limit, is determined only by and The vacuum meta-stability 

condition of is independent of and approximately proportional to [140]. Since 
the loop function is insensitive to for > 100GeV, is almost scaled by 

1/m^^, when the heaviest stau is decoupled. Thus, the excess of is explained by a light 
stau. As found in Fig. 4.8, the upper bound on for larger is more severe than that for 
smaller Such a light stau can be discovered at ILC. 

We discussed that once the deviation from the SM prediction of were observed, the stau 
mass region were determined by the vacuum meta-stability condition. In Sec. 4.3, we will 
study stau searches for ILC based on the stau mass region. 


4.2 Selectron/Smuon 

4.2.1 Overview 

In Sec. 4.1.2 and 4.1.3, we showed that when the muon g — 2 anomaly (2.34) is solved by 
the Bino-smuon contribution (4.5), soft masses of the Bino and the smuon are bounded from 
above by the vacuum meta-stability condition (4.19). The result is summarized in Tab. 4.1. 
In this section, we study experimental status and future prospects to search for such SUSY 
models. The mass bounds depend on the slepton mass spectrum. When the stau is degenerate 
with the selectron and the smuon, a tight constraint is imposed by the vacuum meta-stability 

® The vacuum meta-stability condition determines the upper bound. If thermal transitions are taken into 
account, the constraint could be more severe especially when the stau is light [140]. 

® Although the stau mass region could also be accessed by LHC, future sensitives of the stau searches have not 
been known. In particular, ILC is superior when the stau mass is degenerate with the Bino mass. 
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Mass spectrum 

Smuon 

Vacuum 

LHC/ILC 

LFV/FDM 

— = (Sec. 4.2.2) 

— (Sec. 4.2.3) 

< 330/460 GeV 
< 1.4/1.9TeV 

f 

f or jU 

sf 

sf 


Table 4.1: Summary of searches. The lightest smuon mass is restricted to explain the muon 
g — 2 discrepancy at the lcr/2cj level. The left-right mixing is limited by the vacuum stability 
condition of the stau- or smuon-Higgs potential, depending on the stau masses. The check 
mark shows promising experimental searches. Models of the universal mass spectrum can be 
tested by LHC/ILC, while those of the non-universal spectrum predict large LFV/EDM. 


condition of the staus. In Sec. 4.2.2, it will be discussed that the limit is strong enough for the 
superparticles to be detectable directly in colliders. On the other hand, if the staus are heavier 
or decoupled, the meta-stability bound is relaxed. The superparticle masses can exceed the 
collider sensitivities. In Sec. 4.2.3, it will be argued that such hierarchical mass spectrum can 
be probed by LFV and CPV 


4.2.2 Universal Slepton Mass 

In this section, we study collider searches for the SUSY models with the universal slepton mass 
spectrum. The left- (right-) handed sleptons have a common soft SUSY-breaking mass. 


m 




m 


Mi 


ms 


m 


Sr 


m 


Mr 


m 


Tr- 


(4.24) 


Then, the vacuum meta-stability condition from the stau-Higgs potential restricts slepton 
masses up to 330 — 460 GeV to solve the muon g — 2 anomaly. This is within the reach of 
the LHC or ILC sensitivity. In fact, some of the parameter regions have already been excluded 
by current LHC data, as shown later. 

Collider signatures depend on the LSP In the case of universal soft slepton masses, either 
the lightest neutralino or the lightest stau become LSE^° In the latter case, the stau is likely to 
be long-lived.The (meta-) stable staus leave charged tracks in detectors. Such signatures 
have been studied by ATLAS [149] and CMS [150]. The CMS constraint on the cross section 
of the stau direct production provides the 95% CL exclusion limit, m;^^ > 339 GeV. If this is 
imposed in addition to the vacuum stability condition, all the parameter regions of the long- 
lived stau are excluded in the case of m= < m^o. 

^1 X\ 

In Fig. 4.9, we show the SUSY contributions to the muon g — 2 when the left-right mixing 
is maximized under the two conditions; (i) the vacuum stability and (ii) the neutralino LSI] 

> rtijo. Then, SUSY signature is (opposite-sign same-flavor) di-lepton with large missing 
transverse energy. First, selectrons and smuons are produced by collisions. They subsequently 

The lightest stau is lighter than sneutrinos when fi is large. 

The stau could decay in detectors, for instance, through R-parity violations. The SUSY signatures depend on 
decay channels. Such cases are not discussed here. 
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Figure 4.9: The SUSY contributions to the muon g — 2 as a function of the lightest smuon, 
and the lightest neutralino, m^o. In the orange (yellow) regions, the discrepancy of the 
muon g — 2 (2.34) is explained at the Icr (2cr) level. The left-right mixing is maximized 
under the vacuum stability condition (4.19). The parameters are mj^ = mj^, tan/3 = 40, and 
Mgoft = lOTeV. Further, the condition that the neutralino LSI^ > m^o is imposed. All the 
parameter regions of the long-lived stau < rn^o) are excluded by the ATLAS and CMS 
(long-lived stau search). The region below the green line is excluded by search for direct 
smuon production at ATLAS. 


decay into the lightest neutralino and a partner lepton. This signature was studied by AT¬ 
LAS [118] and CMS [119]. In particular, we show the 95% CL exclusion limit for 
by ATLAS in Fig. 4.9 (and in Fig. 4.4). Here, the region below the green line is excluded. In 
detail, the left-right mixing is negligible in the ATLAS analysis. However, the left- and right- 
handed smuons maximally mix with each other in Fig. 4.9. The total cross section of the 
smuon productions decreases by L0% compared to the ATLAS analysis. Since it is sufficiently 
small, the exclusion limit is considered to be almost the same as the ATLAS result. As a result, 
we find that almost all the Lcr parameter region of the muon g — 2 is already excluded.The 
sensitivity will be improved by the upgrade of the energy and the luminosity. Fig. 4.L0 shows 
the total cross section of the smuon productions at \/s = 8TeV and 14TeV. The cross section 
is estimated at the leading order. For instance, it becomes 1 fb for the lightest smuon mass of 
330 GeV at = 8TeV, which corresponds to the current LHC bound. The same cross section 
is obtained for 450 GeV at L4TeV. Studies for the future sensitivity are required. 


In the ATLAS analysis, is assumed. The selectron production provides almost the same constraint 

as the smuon [118]. Thus, if selectrons are decoupled, the mass bound becomes weaker. However, LFV/CPV 
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Figure 4.10: The total cross section of the smuon productions at LHC with \/s = 8 TeV (dashed) 
and 14TeV (solid). The parameters satisfy mj^ = mj^ = 2Mi. The left-right mixing is maxi¬ 
mized under the vacuum stability condition with tan ^ = 40. 


In the neutralino LSP case, the stau productions have also been studied in LEP and LHC. 
If their mass difference is larger than ISGeV, the stau mass is constrained to be > 81.9 GeV 
by LEP [151]. LHG is still ineffective to search for di-tau events from the direct stau produc¬ 
tions [152]. These limits are so sufficiently weak that Pig. 4.9 does not change even if they 
are imposed. 

In Tab. 4.2, superparticle mass spectra are listed for several points in Pig. 4.9. The lightest 
neutralino mass is close to the Bino mass, since jj, is very large. Given mj = mj^ = mj^, 


constraints are severe, as will be mentioned in Sec. 4.2.3. 



mj 

Ml 




rrin 

TUrj 


^^2 

niyO 


A 

300 

200 

756 

303 

304 

298 

309 

199 

380 

199 

16.1 


300 

200 

699 

303 

304 

299 

308 

209 

375 

199 

14.6 

B 

470 

250 

1680 

472 

472 

465 

479 

329 

581 

250 

10.2 

C 

340 

160 

1138 

343 

343 

336 

350 

199 

442 

160 

18.0 


Table 4.2: Model parameters and mass spectra at several model points in Pig. 4.9. The masses 
are in units of Ge\( and the muon g — 2 is scaled by 10“^°. Here, mj denotes = mj^, and 
tan ^ = 40 is set. 
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Figure 4.11: The production cross section of the lightest smuon at ILC with ^/s = ITeV. 
The beams are assumed to be polarized at 80% (e“) and 30% (e'*'). The left-right mixing is 
maximized under the vacuum stability condition with tan ^ = 40. The Icr (2c7) regions of the 
muon g — 2 are also shown by the orange (yellow) bands. The parameters are = 200 GeV 
and Mgof,. = lOTeV. 


the slepton masses are hierarchical except for those of the selectrons due to a large left-right 
mixing. The lightest stau mass is closest to the neutralino mass among the sleptons. In Fig. 4.9, 
the point A (A) is around the upper side of the contour of the muon g —2. In the vicinity of this 
side, m;? has the closest value to rriyO in the allowed range. If > m^o is imposed, it satisfies 
= m^o to maximize the SUSY contributions to the muon g — 2 around the upper side of 
the contour. Above it, the regions are excluded by the long-lived stau search, or the muon 
g — 2 becomes too small. On the other hand, if the mass difference between the stau and the 
neutralino is assumed to be larger than 5M, they satisfy = m^o -I- 5M in the vicinity of the 
upper side. At the point A, 5M = 10 GeV is imposed, and is larger than m^o by 10 GeV. 
In contrast, the points B and G are away from this region. B (G) is close to the maximal end 
point of the lightest smuon mass which explains the muon g — 2 at the 2cr (Icr) level. Here, 
the left-right mixing is determined by the vacuum meta-stability condition of the stau. 

Linear colliders of (ILG) are very useful (see e.g.. Ref. [148,153]). They can provide 
rich informations of the models. Moreover, they are superior to LHG when the slepton masses 
are close to that of the lightest neutralino [154]. In the orange/yellow regions of Fig. 4.9, 
selectrons, smuons and staus can be produced at the linear colliders. In Fig. 4.11, we show 
production cross sections of the lightest smuon for the left- and right-handed smuon masses. 
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Here, the production cross section of the lightest smuon are [155] 
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(4.25) 


and it is assumed that the collision energy is \/s = 1 TeV at ILC, and the beams are polarized at 
80% for the electron and 30% for the positron. The parameters are defined in Appendix C. The 
left-right mixing is maximized under the vacuum stability condition. The other parameters are 
Ml — 200GeV and tan/3 = 40. It is found that the cross section is larger than 20 (l)fb“^ for 
the Icr (2cr) parameter region of the muon g — 2. It is expected that the smuons can be 
discovered in almost all the parameter region that are kinematically allowed (cf. [148]). 

It is possible to measure masses of the smuon and the neutralino from event distributions 
and the cross section [124,156]. Also, the chirality structure of the smuons could be deter¬ 
mined by the beam polarization. We assume {P^-,P^+) — (+0.8,—0.3) in the left panel of 
Fig. 4.11, and (—0.8, +0.3) in the right panel. The cross section of a chiral smuon is sensitive 
to the polarization, because the productions proceed by the s-channel y/Z exchanges. On the 
other hand, when = rriji^, the cross section is insensitive to the beam polarization. Since 
the left-right mixing is large, the left- and right-handed smuons are maximally mixed with 
each other. Both of them can be produced, and the mass difference is expected to be mea¬ 
sured, for instance, by measuring threshold productions. At the model points in Tab. 4.2, the 
difference is > 10 GeV, which is much larger than the uncertainty of the mass measurement 
at ILG, ^(10 — 100) MeV [153]. It is emphasized that the smuon productions are clean and 
direct signatures of the SUSY contributions to the muon g — 2. 

The selectron productions proceed not only by the s-channel y/Z exchanges, but also 
by the t-channel Bino exchange. The latter contribution can enhance the cross section and 
provide individual information in addition to mass measurements of the selectron and the 
neutralino. The productions, e^e~ —> e^e~ and e+e“ —> proceed by the t-channel. In 

Fig. 4.12, we show their total cross sections by the black solid lines (For the formulae of 
the cross section, see Appendix G.). Here, the parameters are and the beam 

polarization is 80% for the electron and 30% for the positron. It is found that the cross 
sections are ^(10) fb in the muon g — 2 parameter regions. They are larger than those of 
the smuon in Fig. 4.11. Importantly, this channel is useful to measure the Bino-electron- 
selectron couplings [124,125]. They are deviated from the U(l)y gauge coupling constant by 
decoupling heavy superparticles, as discussed in Sec. 4.1.2. In Fig. 4.12, the cross sections are 
estimated with (without) dg^ and Sg^ (see Eqs. (4.11) and (4.12)). The results are shown by 
the solid (dashed) lines. It is found that the corrections enhance the cross section by 8 — 10% 
for Mgusy = 10 TeV. This is measurable at ILG [124,125]. 

The cross section is calculated at the leading order in Fig. 4.11. It can be enhanced by ~ 10% near the mass 
threshold [155], 
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Figure 4.12: The production cross section of e+e“ and at ILC with ^/s = ITeV. The 
parameters are same as Fig. 4.11, but is varied with = ntg. The corrections to 

the Bino coupling due to decoupling of heavy superparticles are included (discarded) in the 
black solid (dashed) lines. 


Cross sections of e^e“ —» and e^e~ —» e’^e“ also include the t-channel contribution 

of the Bino. In particular, the former cross section is enhanced well compared to those solely 
by the s-channel y/Z exchanges. It can be > 100 fb when is relatively small in the muon 
g — 2 parameter regions. However, in the latter process, the t-channel contribution interferes 
destructively with the s-channel contribution. On the other hand, the cross section of e~ —> 
e^e“ differs by 5 — 6% between the cases with and without and 5g^. This is smaller than 
the above channels. 

Let us comment on the dark matter in the present model. The lightest neutralino is a 
candidate of the dark matter. The neutralino relic abundance becomes consistent with the 
measured cold dark matter abundance [157,158] by the stau co-annihilation. The mass dif¬ 
ference between the stau and the neutralino is required to be 5 — 10 GeV. This corresponds 
to the model points in the vicinity of the upper side of the muon g — 2 contours in Fig. 4.9, 
including the points A or A’. This region has not been excluded by the studies of the di-lepton 
signature at LHC [118,119]. Masses of the selectrons and the smuons are too close to that of 
the neutralino for LHC. ILC is superior to study the region [148,153]. Selectrons and smuons 
as well as staus can be produced, as discussed above. 

Finally, we comment on the Higgs coupling to di-photon. As seen in 4.1.4, when the stau 
is light and the left-right mixing is large, the coupling can be sizably enhanced. Comparing 
with the muon g — 2 parameter regions in Fig. 4.9, the stau contribution 5 becomes 5 — 20% 
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Figure 4.13: The lower bound on the stau mass to satisfy the vacuum stability bound of the 
stau-Higgs potential. In the orange (yellow) region, the SUSY contribution to the muon g — 2 
can explain the discrepancy (2.34) at the Icr (2cr) level. The parameters are = M^, 

tan^ = 40 and = 30TeV. 


(2 — 5%) in Icr (2cr) region. The latter is comparable with or larger than the future sensitivity 
to using joint analysis of HL-LHC and ILC. If the deviation is observed by future collider 

experiments and can be explained by the stau contribution, we will get more strong evidence 
for the SUSY. The stau searches are discussed in Sec. 4.3. 


4.2.3 Non-Universal Slepton Mass 

In this section, we discuss the case without slepton mass universality. The slepton soft SUSY- 
breaking masses satisfy 


(4.26) 

As discussed in Sec. 4.1.3, the vacuum meta-stability bound of the stau-Higgs potential is 
relaxed when the staus are heavy, and the smuon masses are allowed to be larger in order 
to solve the muon g — 2 anomaly. In other words, the stau masses are required to be large 
relative to the smuons, when the smuons are heavy. In Fig. 4.13, we show the lower bound 
on the stau mass. In the orange (yellow) regions, the muon g — 2 anomaly is solved at the Icr 
(2cr) level, while the vacuum meta-stability constraint is avoided. The ratios of the stau and 
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smuon masses are defined as 
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tJ-R 


The parameters are tan^ = 40 and = SOTeV. The vertical 

axis is R = R^ — Rj^. If the staus are decoupled, the left-right mixing is bounded either by 
Eq. (4.21) or the vacuum meta-stability bound of the smuon-Higgs potential, as discussed in 
Sec. 4.1.3. The smuon masses can be 1.4TeV (1.9TeV) for the la (2a) level of the muon 
g — 2. Even though it becomes difficult to produce such sleptons at LHC or ILC, non-universal 
slepton masses generically cause problems of too large lepton flavor violations (LEV) and CP 
violations (CPV). 

When the muon g — 2 anomaly is solved by the SUSY contributions, LEV and lepton EDM 
are generically sizable. In this section, the following setup is considered. In the model basis, 
the slepton mass matrices are diagonal among the flavors. 
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(4.28) 

(4.29) 


Eurthermore, mg^ = and mg^ = are imposed. Otherwise, LEV becomes too large, 
as shown later. The Yukawa matrices are generally non-diagonal in the model basis. The 
mass eigenstate basis of the charged leptons is obtained by the left- and right-handed unitary 
matrices, U^, Ur, as Eq. (3.19). In this dissertation, we assume that the unitary matrices are 
represented as 


(4.30) 


Non-vanishing mixings, (5^^ r);j, induce LEV and EDM, as long as the slepton masses are non- 
universal.^"^ On the other hand, we assume that quark ECNC and CPV are suppressed by heavy 
colored superparicles in this dissertation. 

The magnetic dipole contributions to LEV are represented in Eq. (2.10), as mentioned in 
Sec. 2.1.2. In our setup, the Wilson coefficients. At. and At^., are dominated by the Bino- 
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slepton contributions, similarly to the muon g 
are estimated as [2,17] 


2. In the mass insertion approximation, they 
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Similar setup has been studied in Ref. [159], where squarks of the first two generations are light. 
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The two loop factor is found in Eq. (4.9). The loop function F^j, is defined as 
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mt 
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■ a,b 


2 2 

m~ m~ 

ij.n il 


-fl 


(■Rb 


N 


W’ |MJ 


(4.33) 


The muon FCNCs are most sensitive to the non-universal slepton mass and the non-diagonal 
Yukawa matrix. In this case, the coefficients is important. The unitarity matrices in A ^2 
are expanded as 




a,b=l,2,3 


m 




1 +A, 


m, ^ N 

I * (^l)i3(^r)23 (^1,2 ~ -^1,3 ~ ^ 3,2 + -^3,3) 3 


■(^l)i2 (^1,2 ^ 2 , 2 ) 


(4.34) 


1 +A 


at the leading order of (5^);^, (5^)ij and Here and hereafter, we set = 0, for 

simplicity. Similarly, A ^2 is obtained by replacing L R. In the last term, 2 is dominant 
when the staus are heavy, i.e., F 12 » F 1 , 3 ,F 3 , 2 ,P 3,3 for ^ ^ ^^r- 

On the other hand, the right-hand side vanishes when the slepton masses are universal, as 
expected from the super GIM mechanism. In particular, when the sleptons are degenerate in 
the first two generations, the first term becomes zero due to 2 = ^ 2 , 2 - Otherwise, the muon 
LEV is induced at the order of ( 5 ^) 12 . 

The decay rate of /t —> ey is evaluated by Eq. (2.39). It is compared with the SUSY 
contribution to the muon g — 2. In the non-universal slepton mass spectrum, it is represented 
byA^. andA^. as 

j i] i] 


a„(SUSY) = m=Re[A‘,+A*,] 


a 


(1 -I- —m„MijUtan B 

471 


m 


V- 


1 +A 


F2,2 + ^ 




This is same as Eq. (4.5) up to a correction k, which is represented as 




K = 


1 +A, 


-Re [(5i)23(5^)23] (^ 2,2 F 23 R3.2 + R3,3)+' 


(4.35) 


(4.36) 


Here, the omitted terms are suppressed by higher order of (5^);^, (5^^);^ or m^/m^. If the 
slepton mass matrices are universal, k vanishes. Eor the mass spectrum (4.26) with » 1, 

the ratio is 


Br(/x ^ ey) , 1 I c r |2 f 1 +A^ y 

a^(SUSY)2 16 I '' "'I l,m^l + Aj ’ 


(4.37) 


where Tj^j is the total decay rate of muon, and the mixing is defined as 

















86 


CHAPTER 4. PROSPECTS EOR SLEPTON SEARCHES 


CO 

*0 

CO 

tH 

'a 



10"^ 


—1- 1 - 1 —• 1 • 1 -; 

10"® 

l\ \ 

I \ 

a^(SUSY) = 10"^ : 

10"® 

1 \ 

4 \ 

4e(YbF) , 


\ 

' \ \ 

• \ 

\ - 

ThO - 

10“® 

\ 

. \ \ 

• \ N 

■ \ 

Pr, ThO 

10“® 


YbF,WN - 

10-10 


.1 


1 10 


R = mf/mfx 


Figure 4.14: Contours of Br(;U —> ej) (left) and (right) for a^(SUSY) = 10“^ with = 
SOTeV. In the left panel, the upper black solid line is the current bound by MEG, Br(/t —> 
ey) < 5.7 • 10“^^ at 90% CL [50] . The lower one is the sensitivity of the MEG upgrade, Br(jU —» 
ey) — 6 • 10“^^ [53]. The blue dashed line is the Mu3e sensitivity, Br(jU —> eee) = 10“^^ [55]. 
The jJL — e conversion is expected to probe down to R^^ = 3 • 10“^^ (upper red dotted) by 
COMET/Mu2e [57,58] and 2 • 10“^^ (lower red dotted) by PRISM/PRIME [59]. In the right 
panel, the black solid line is the bounds by using YbE and ThO, |dj < 1.05 • 10“^^ecm and 
8.7X 10“^^ecm at 90% CL [61,62]. The sensitivity is planned to be improved: |dj = 10“^^ ecm 
by Er or ThO [64-66], and 10“^° ecm by YbE or WN [63,67]. Here, rtii^ = mi^ is assumed. 


It is independent of superparticle mass spectra except through and A.J.. Thus, when the 
muon g — 2 discrepancy (2.34) is explained by the SUSY contributions, jji ey is induced 
by the non-universal slepton mass sizably. It is important that the decay is not suppressed 
by heavy slepton masses, for instance, = 1.4TeV or 1.9TeV in Tab. 4.1, for given SUSY 
contributions to the muon g — 2. 

In the left panel of Eig. 4.14, contours of Br(ju ^ ey) are shown. Here, is 

assumed. The SUSY contributions to the muon g — 2 is fixed to be a^(SUSY) = 1 x 10“^ 
with Mgoft = 30 TeV. The small corrections are taken at = 400 GeV and tan ^ — 40 

as a reference, though the result is almost independent of them. In the figure, the contours 
correspond to the current limit and future sensitivities of experiments, 

• the current limit of the MEG experiment, Br(;i —> ey) < 5.7 x 10“^^ at 90% CL [50] 
(upper black solid line in the figure). 

• the sensitivity of the MEG upgrade, Br(jU —> ey) = 6 x 10“^"^ [53] (lower black solid 
line). 
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• the sensitivity of the Mu3e experiment, Br(jH —> eee) = 10 at Phase II [55] (blue 
dashed line). 

• the sensitivity of the COMET experiment, = 3 x 10“^^ at Phase II [57] (upper red 
dotted line). The Mu2e experiment has a similar sensitivity [58]. 

• the proposal of the PRISM/PRIME project, R^^ = 2 x 10“^^ [59] (lower red dotted line). 

Note that jjl —> eee and jjl — e conversion experiments have better sensitivity in future than 
those of jU —» ey. In particular, the latter experiment has low (accidental) backgrounds. On 
the other hand, the current constraint and future sensitivities of the tau LEV’s are weaker than 
those of the muon, though they are also induced by the non-universal slepton masses with 
finite ( 5 J 23 and ( 5 ^) 23 . 

Erom the figure, it is found that the LEV decay rate increases rapidly as the staus become 
heavier than the smuons, R > 1. When the staus are decoupled, R » 1, the mixing 1 5 ^ 3523 1 
in Eq. (4.38) is limited to be smaller than 3 x 10“^ by MEG for a^(SUSY) = 1 x 10“^. Eor 
instance, when the smuon and selectron masses are larger than 1 TeV, the stau are required to 
be heavier than 7TeV (R > 7) to avoid the vacuum stability constraint, according to Eig. 4.13. 
If the lepton Yukawa matrix is related to the quark sector, e.g., by the GUT relation, it is 
naively expected to be |5i3 523l ~ KbKb ~ 10“^^. This already exceeds the above limit. Thus, 
non-universal slepton mass spectra are tightly constrained by LEV In future, if sleptons are 
neither discovered at LHG nor ILG, the model is expected to be probed by LEV Otherwise, 
the SM Yukawa matrices are tightly limited in the model basis, or when flavor models are 
constructed. 

The flavor off-diagonal components of the Yukawa matrices are sources of the GP viola¬ 
tions, when the slepton mass matrices are non-universal. Similarly to the correction k of the 
muon g — 2 in Eq. (4.35), the electron EDM is induced as 



= (1 + 52'°°P)—Mi/rtan^ 


871 




Here, we assume jS) = 0. The omitted terms are suppressed by orders of 

(^R)ij OP In tfio last term, 3 » F3 3, F3 3, F3 3 is obtained when the staus are heavy. On 

the other hand, the right hand side vanishes when the slepton masses are universal, because 
the complex phases can be rotated away. Gomparing Eq. (4.39) with Eq. (4.35), we obtain 


dje 1 



(4.40) 


a^(SUSY) 2m^ 
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for the mass spectrum (4.26) with » 1. It is independent of superparticle mass spectra 

except through and A^. Thus, if the muon g — 2 anomaly is solved by the SUSY contribu¬ 
tions, EDM becomes sizable by the non-universal slepton mass. Similarly to Lpy this is valid 
even for large selectron and smuon masses. 

In the right panel of Fig. 4.14, contours of the electron EDM are shown. Here, the SUSY 
contributions to the muon g — 2 is fixed to be a^(SUSY) = 1 x 10“^ with = SOTeV. The 
result is almost independent of superparticle mass spectra except for small corrections, A^ 
and A^. The following data are used, 

• the limit with the YbF molecule, \dj < 1.05 x 10“^^ecm at 90% CL [62] (upper solid 
line in the figure). 

• the current limit with the ThO molecule, \d^\ < 8.7 x 10“^^ecm at 90% CL [61] (lower 
solid line in the figure). 

• the sensitivity with the Fr atom, \dj = 10“^^ecm [64] (dashed line). The experiment 
with the ThO molecule could have a similar sensitivity by accumulating data, \dj = 
1 X 10“^®/ (day) ecm [65,66]. 

• the sensitivity with the YbF molecule, IdJ = 10“^° ecm [63] (dotted line). The experi¬ 
ment with the WN ion can probe down to, |dj = 10“^° ecm/day, where the systematic 
limit is at the level of 10“^^ ecm [67]. 

From the figure, it is found that EDM is sensitive to Im[(5^)^3(5i)]3] when the staus are heav¬ 
ier than the selectrons. The current experimental limit puts a constraint, Im[(5^ )i3(^l)i3] < 
6 X 10“® for R:> 1 and a^(SUSY) = 1 x 10“^. The sensitivity will be improved very well. The 
mixing will be able to be probed at the level of 10“^°. Thus, if sleptons are neither discovered 
at LHC nor ILC in future, the model can be sensitively probed by EDM as well as LFV If no 
signal will be observed, the CP violating phase must be suppressed very tightly in order to 
explain the muon g — 2 anomaly by the SUSY contributions. 


4.3 Stau 

So far, we studied the prospects for the selection/smuon searches. In this section, we discuss 
the stau searches. If the staus have masses of ^(100) GeV and large left-right mixing, the 
Higgs coupling to the di-photon is deviated from the SM prediction. Such a large left-right 
mixing is constrained by the vacuum meta-stability condition of the stau-Higgs potential, as 
mentioned in Sec 4.1.3. Then, we discussed that once the deviation from the SM prediction of 
were observed, the mass region for the staus were determined by the condition in Sec. 4.1.4. 

Once the stau is discovered at ILC, its properties including the mass are determined. Par¬ 
ticularly, it is important to measure the mixing angle of the stau 6^. When sin 20^ is sizable, 
the angle can be measured at ILC [124,160-162]. As observed in Fig. 4.8, it is likely to be 
sizable to enhance k^. In particular, if sin20;f is large enough to be measurable, the heaviest 
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Table 4.3: Model parameters at our sample point. In addition, tan ^ = 5 and = 0 are set, 
though the results are almost independent of them. 


Parameters 


^^2 

sin 20^ 

rriyO 

5k.^ 

Values 

100 GeV 

230 GeV 

0.92 

90 GeV 

3.6% 


stau is likely to be light. Thus, it may be possible to discover the heaviest stau and measure 
its mass at ILC. Then, the stau contribution to jc^, can be reconstructed by using the measured 
masses and mixing angle. This provides a direct test whether the stau contribution is the ori¬ 
gin of the deviation of /c^. On the other hand, the heaviest stau is not always discovered at the 
early stage of ILC, even if the stau mixing angle is measured. If 6^^ as well as is measured, 
may be estimated in order to explain the excess of fc^. In this section, we will study the 
reconstruction of the stau contribution to fc^. The mass of the heaviest stau and theoretical 
uncertainties will also be discussed. 


4.3.1 Reconstruction 

If both of Ti and T 2 ^re measured, the stau contribution to can be reconstructed. The 
contribution is determined by the parameters in Eq. (4.22). In this subsection, we assume the 
situation that both of Tj and T 2 are observed by ILC. Then, we discuss how and how accurately 
the model parameters are measured at ILC, and consequently the stau contribution to is 
reconstructed. 

Let us first specify a model point to quantitatively study the accuracies. In Tab. 4.3, the 
stau masses, the stau mixing angle, and the Bino (-like neutralize) mass are shown. The 
point is not so far away from the SPSla’ benchmark point [163], where ILC measurements 
have been studied (see e.g.. Ref. [160]). The stau mixing angle is chosen to enhance the 
Higgs coupling as = 3.6%. The staus masses are within the kinematical reach of ILC 
at ^/s = 500 GeV. The point is consistent with the vacuum meta-stability condition and the 
current bounds from LHC and LEE The most tight bound on the stau mass has been obtained 
at LEE as > 81.9GeV at 95% GL [151]. On the other hand, LHG constraints are still 
weak [152]. The other superparticles are assumed to be decoupled for simplicity. In particular, 
tan p = 5 and A^ = 0 are chosen, where the Higgsino masses are about 2.2 TeV. 

In order to reconstruct the stau contribution to it is required to measure the stau 
masses and the mixing angle. At ILC, staus are produced in e~^e~ collisions and decay into 
the tau and the Bino. The stau masses are measured by studying the endpoints of the tau 
jets. In Ref. [160], the mass measurement has been studied in detail at SESla’. It is ar¬ 
gued that the mass can be measured at the accuracy of about 0.1 GeV (6GeV) for Tj (T 2 ). 
Here, ^/s = 500GeV and .Sf = 500fb“^ are assumed for ILG. The mass resolution may be 
improved by scanning the threshold productions [153,164]. The accuracy could be ~ IGeV 
for 771;^^ = 206 GeV. Since the model parameters of our sample point are not identical to those 
of SESla’, the mass resolutions may be different from those estimated at SESla’. Eor instance. 
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Figure 4.15: Contours of 5 sin20;f/sin20^ determined by the measurement of the production 
cross section of a pair of (left) and that of and T 2 (right). Uncertainties from the mass 
resolutions are not taken into account. 

the production cross section of staus becomes different, while the SUSY background is negli¬ 
gible in our sample point. Profile of the tau jets depends on the masses of the staus and the 
Bino. In this dissertation, instead of analyzing the Monte Carlo simulation, we simply adopt 
the mass resolution, 

~ 0.1 GeV, ~ 6GeV. (4.41) 

Next, let us discuss the measurement of the stau mixing angle, 9;^. Several methods have 
been studied for ILG. For instance, the polarization of the tau which is generated at the stau 
decay has been studied in Ref. [124,160,161]. The angle can also be extracted from the 
production cross section of a pair of the lightest stau [161]. Note that accuracies of these 
angle measurement depend on the model point, i.e., the input value of 0^. 

In order to study the accuracy of the stau mixing angle at our sample point, let us investi¬ 
gate the production cross section of the lightest stau by following the procedure in Ref. [162]. 
The production cross section of the lightest stau is given by [161] (cf. Appendix G) 

^ __^ 8na^ 3 r , 

CT(e+e ^TiTi) = A2 Cii . 4 (P-+L + P+-R ) 

3s L sm 20^v 

The resolutions estimated in Ref. [160] depend on the uncertainty of the measured Bino mass. The Bino 
mass can be measured very precisely at ILC by the productions of selectrons or smuons [148], though they are 
irrelevant for and the vacuum meta-stability condition. 
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1 

+ 77 (^-+ + ^+-) + ^11 


2 sin^ 20 ^ 1 / 


(P_+L + P+_il) 


(4.42) 


at the tree level, where the parameters are A = 1 — 4m~ /s, = s/(s — mp, and = 

[(L+P) + (L —]?)cos20^]/2 with L — —1/2 + sin^ 9^^ and R = sin^ 9^^. The beam polarizations 
are parameterized as P^± = (1 TPe-)(l ^Pe+)- In the bracket, the first and second terms come 
from the s-channel exchange of the Z boson and the photon, respectively. The last term is 
induced by the interference of them. The dependence on the stau mixing angle originates in 
the Z boson contribution. 

Since Eq. (4.42) is a function of the stau mass and mixing angle, 9:^ is determined by 
measuring the cross section and the stau mass. If the errors are summed in quadrature, the 
accuracy of the stau mixing angle is estimated as 


(Asin20;f)^ 


d sin20;f h ^ 

5cj(ti) J 


(Acr(Ti))V 


d sin 2 0s 


dm 






(4.43) 


where cr(Ti) = (j{e~^e~ —> In Fig. 4.15, we show contours of the uncertainty of the stau 

mixing angle, 5 sin 2 0;^/ sin 20;^. In the left panel, the angle is determined from the production 
cross section of the lightest stau. The accuracy is sensitive to the input value of sin20;f and 
5cr(Ti)/cj(T;^). In contrast, the uncertainty from the mass resolution of in Eq. (4.41) is neg¬ 
ligible. The accuracy of sin20;f becomes better for larger sin20;f. If the stau contributes to 
sizably sin 20^ is likely to be large, as observed in Fig. 4.8. Thus, the mixing angle is expected 
to be measured well. At the sample point, where sin20;f = 0.92, 5 sin20;f/sin20;f is estimated 
to be better than 10%, if the cross section is measured as precisely as 5cr(Ti)/cr(T^) < 10%. 
At ILC, it is argued that the production cross section can be measured at the accuracy of about 
3%, according to the analysis in Ref. [160] at SPSla’. If 5c7(Ti)/cr(T^) ~ 3% is applied to our 
sample point, the accuracy is estimated to be 


Asin20;jr/sin20;f ~ 2%. 


(4.44) 


From Eqs. (4.41) and (4.44), the accuracy of the reconstruction of the stau contribution to 
Ky is estimated. If the errors are summed in quadrature, the uncertainty is estimated as 



(A sin 20^) ^ (4.45) 


Substituting Eqs. (4.41) and (4.44) for Eq. (4.45), the uncertainty is obtained as 


Ak^ ~ 0.5%, (4.46) 

at the sample point, where = 3.6%. Note that the uncertainty of the measurement of 
fc^ is 1-2% from HL-LHC and ILC, as mentioned above. Since the reconstruction error is 
comparable to or smaller than that of the measured it is possible to check whether the stau 
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is the origin of the excess of the Higgs coupling k^. It is emphasized that this is a direct test of 
the Stan contribution to k^. 

In Eq. (4.46), the error is dominated by the uncertainties of the heaviest stau mass and 
the stau mixing angle. The former may be reduced by scanning the threshold of the stau 
productions, as mentioned above. For instance, if we adopt ~ 1 GeV as implied in 

Ref. [153,164], the error becomes Afc^ ~ 0.3%. On the other hand, the latter uncertainty 
may be improved by studying the production cross section of and T 2 [160]. The cross 
section is given by [161] (or see Appendix C) 


cr(e+e 


Sna^ 3 

-A2 

3s 


2 X c 


single 


-(P_+L^ + R+_R^) 


w 


(4.47) 


where A = [1 - (m^^ + m;fjVs][l - (m^^ - and c ^2 = C 21 = (L - R)sin29:^/2. 

Since —> T 1 T 2 proceeds by the s-channel exchange of the Z-boson, its cross section is 
proportional to sin^26;f. Thus, it is very sensitive to the stau mixing angle, and further, the 
accuracy is independent of the model point, once the error of the production cross section is 
given. In the right panel of Fig. 4.15, we show contours of 5 sin20;f/ sin20;jr, where, the angle 
is extracted from a(e~^e~ —» T 1 T 2 ). It is found that the accuracy is independent of the input 
sin20;f. Here, uncertainties from the mass resolution are neglected. In particular, if the mass 
resolution of T 2 is large, the accuracy of the mixing angle becomes degraded. Unfortunately, 
the accuracy of the measurement of C7(e'''e“ —> t^T 2 ) has not been analyzed for ILC. Since 
sin20;f is likely to be large to enhance k^, the cross section can be sizable. At the sample 
point, it is estimated to be about 6fb for -/s = 500GeV with {P^_,P^_^) = (—0.8,0.3). It is 
necessary to study this production process in future. 

Let us comment on the tan/3 dependence. At the sample point, we set to tan^ = 5. 
Although the stau contribution to includes tan^, once the stau masses and mixing angle 
are measured, the reconstruction of the Higgs coupling is almost insensitive to it. This is 
because the stau left-right mixing parameter m~ is determined by the measured masses and 
mixing angle through Eq. (4.23). One can check that, even if tan^ is varied, the accuracy 
(4.46) is almost unchanged. 

Finally, we comment on the muon g — 2. In this section, we discussed our sample point, 
where rri:^^ = 100GeV, = 230GeV, sin20;f = 0.92, m^o = 90GeV, and tan^ = 5, respec¬ 
tively. The soft SUSY breaking masses can be estimated as = 145 GeV, = 291 GeV,^^ 
and the SUSY contribution to the muon g — 2 is a^(SUSY) ~ 1.8 x 10“^, if the sleptons masses 
are universal among the flavors. Here, we set to = 30TeV and assume that the left-right 
mixing of the smuon is = (m^/m^) x It is found that our sample point is explained 
the muon g—2 discrepancy at Icj level and is not excluded by dilepton searche [149,150] . The 
mass difference between and m^o are ~ 50 GeV, which is much larger than the uncertainty 
of the mass measurement by smuon productions at ILC. Therefore, the smuon productions are 
expected to provide direct signatures in this scenario. 


In this case, we consider 0 < 0^ < f region. 
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4.3.2 Discussions 

Let us discuss miscellaneous prospects for the stau searches and theoretical uncertainties 
which have not been mentioned so far. 

First of all, let us consider the situation when the heaviest stau is not discovered at the 
early stage of ILC, i.e., at ^/s = 500 GeV. As found in Sec. 4.1.4, if the excess of is measured 
at this stage, the lightest stau is already discovered. Then, it is possible to determine the stau 
mixing angle by measuring the production cross section of the lightest stau, as long as sin 26^ 
is sizable as observed in Fig. 4.15. From the measurements of m;^^, 6^ and jc^,, the mass of the 
heaviest stau is determined. The predicted mass could be tested at the next stage of ILC, 
e.g., ^/s = ITeV. 

In order to demonstrate the procedure, let us consider a model point with = 150 GeV, 
= 400 GeV and sin20;f = 0.91. At the point, the deviation from the SM prediction of the 
Higgs coupling is = 5.6%. At the early stage of ILC, it is expected that is determined 
with the uncertainty A/c^ ~ 2%, and the lightest stau is measured with Am;^^ ~ 0.1 GeV and 
5cr(Ti)/cr(Ti) 3%. The stau mixing angle is extracted from the cross section as 

A sin 26^/sin 20^ ~ 2.5%, (4.48) 


as observed in Fig. 4.15. Since 5 k:^ is a function of rriT^^ and 0^, the uncertainty of the 
heaviest stau is estimated as 




dm 


-^2 


dm^ 


(Am^j" + 


dm^f 


dm 


d sin 20s 


(Asin20^)"+I I (AkJ ,(4.49) 


thus the mass of the heaviest stau is determined with the accuracy 


Am^^~53GeV, (4.50) 

where, the largest uncertainty comes from the measurement of the Higgs coupling. If the 
error is reduced to be Afc^ ~ 1% due to reductions of the HL-LHC systematic uncertainties 
(see Ref. [74]), Am;^^ ~ 26GeV is achieved. Such a prediction can be checked at ILC with 
^/s = 1 TeV. This result would be helpful for choosing the beam energy to search for the 
heaviest stau at ILC. Once T 2 is discovered, the stau contribution to the Higgs coupling can be 
reconstructed as Sec. 4.3.1. 

When the sleptons have universal mass spectrum, the left- and right-handed slepton masses 
can be estimated by mj = 345 GeV, mj — 250 GeV, respectively.^^ If the lightest neutralino 
are much lighter than the smuons, they are expected to be probed at LHC with ^/s = 13- 
14TeV. The SUSY contributions to the muon g — 2 is calculated as a^(SUSY) ~ 2.0 x 10“^, 
and is explained the discrepancy of the muon g — 2 at Icj level, where the parameters are 
taken to = 140 GeV, and = 30 TeV. Combining the information of the Higgs coupling 
with that of the muon g — 2, we may detect effects of sleptons more definitely. 


We consider ^ < 6^ < ^ region. 
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Figure 4.16: Chargino contribution to the Higgs coupling The left panel shows contours 
of the chargino contribution to where the Wino mass M 2 is equal to the Higgsino 

mass jJL. In the right panel, the chargino contribution 5k^(x^) is displayed for various tan^, 
where is fixed to be 250 GeV. 


The uncertainty of the prediction of the heaviest stau mass depends on the model point, 
especially the stau mixing angle. If is larger, sin20;f is likely to be smaller, as expected from 
Fig. 4.8. The measurement of the stau mixing angle, then, suffers from a larger uncertainty, 
and it becomes difficult to determine the mass of the heaviest stau. 

Next, let us mention extra contributions to the Higgs coupling from other SUSY particles. 
So far, they are suppressed because those particles are supposed to be heavy. However, if 
the chargino, the stop or the sbottom is light, its contribution can be sizable [95,144]. These 
particles are searched for effectively at (HL-) LHC (see e.g.. Ref. [146,147]).^® If none of them 
is discovered, their masses are bounded from below, and upper limits on their contributions 
to Ky are derived.These extra contributions should be taken into account as a theoretical 
(systematic) uncertainty in the analysis of 5 k^. 

We discuss the chargino contribution. Since the mass bounds on the charginos are still 
weak, it can be as large as the theoretical uncertainty. At one loop level, the formula was 
given in (3.65) of Sec. 3.5. When the charginos are heavier than the electroweak gauge 


On the other hand, it is possible to determine tan [3 by studying decays of staus, neutralinos or charginos at 
ILC, if the Higgsinos are light [153,161]. 

If extra SUSY particles such as charginos are discovered, their contributions to k^. may be reconstructed. 
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bosons, it is approximated as [144] 




4 g^vsin2/3 
3 M 2 ^ - sin2l3 


(4.51) 


Since sin2^ = 2tan^/(l +tan^/3), it is found that it is suppressed in large tan^, thus is 
relevant when tan ^ = ^(1). 

The left panel of Fig. 4.16 shows contours of Sk.^^x'^) as a function of tan^ and the 
lightest chargino mass my±. Here, we take M, = jU, for simplicity. The contribution decrease 

A 1 

as my± or tan B increases. If the charginos are constrained to be heavier than 600 GeV (1 TeV), 

A 1 

these contribution to is estimated to be smaller than 0.5% (0.2%) for tan^ > 2. This is 
considered to be a theoretical uncertainty. In addition, if either of the Wino or the Higgsino is 
decoupled, is suppressed. Such a feature is observed in the right panel of Fig. 4.16, 

where is fixed to be 250 GeV for various tan^. 

The stop and the sbottom can also contribute to fc^ sizably [95]. In addition, they simul¬ 
taneously modify the Higgs coupling to di-gluon. The coupling is expected to be measured 
precisely at ILC at (sub) percent level [74,75]. Hence, if deviations are discovered in as 
well as Ky, it is interesting to study the contributions of stop or sbottom. 
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Chapter 5 
Conclusion 


Supersymmetry is one of the most motivated candidates of new physics. If the muon g — 2 
anomaly is solved by SUSY, some of superparticles are relatively light. In this dissertation, we 
concentrated on the "minimal" SUSY model in which only the Bino and the slepton are light, 
while the other superparticles are decoupled. In this model, the muon g — 2 discrepancy is 
explained by the Bino-smuon diagram. It is enhanced by large left-right mixing m~ . We have 

^LR 

investigated various phenomena in the minimal model. 

The analyses were classified by the slepton mass spectrum. When the mass spectrum 
is universal among the flavors, we found that the vacuum meta-stability of the stau-Higgs 
potential restricts the smuon masses tightly. They are predicted to be within 330 (460) GeV 
at the Icr (2cr) level of the muon g — 2. It was shown that part of the parameter region is 
already excluded by LHC, and argued that such sleptons are expected to be probed at LHC or 
ILC in future. 

If the staus are (much) heavier than the smuons, the vacuum stability bound of the staus 
is relaxed. In this mass spectrum, the smuon masses are limited by the vacuum mata-stability 
condition of the smuon-Higgs potential. It was found that they are less than 1.4 (1.9) TeV 
at the lc7 (2cr) level of the muon g aLS 2. Such slepton masses exceed the LHC/ILC reach. 
Instead, the non-universal slepton mass spectrum generically predicts too large LFV and EDM. 
They originate in the non-diagonal SM Yukawa matrices in the model basis. Although the 
prediction depends on the flavor models, LFV and/or EDM is very likely to be sizable in wide 
models. If no sleptons are discovered at LHC or ILC, the model is expected to be probed by 
LFV and EDM sensitively. 

In the universal mass spectrum case, the stau constribution to the Higgs coupling to di¬ 
photon becomes large. The coupling will be measured at the percent levels by the joint 
analysis of HL-LHC and ILC. In this dissertation, we considered a situation that an excess of 
is measured in HL-LHC and ILC. First, we studied the stau mass region based on the vacuum 
meta-stability of the stau-Higgs potential. Consequently, we showed that the lightest stau 
is predicted to be lighter than about 200 GeV, if the excess is measured to be larger than 4% 
in the early stage of -^s = 500 GeV ILC. Further, it was found that is limited to be less than 
290-460 GeV, if the excess is measured to be 1-2% by accumulating luminosity at -v/s = 1 TeV 
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ILC. Such stau is within the kinematical reach of ILC. Therefore, we concluded that the stau 
contribution to fc^, can be probed by discovering the stau, if the excess is measured in the 
future experiments, and if it originates in the stau contribution. 

Once the stau is discovered at ILC, its properties are determined precisely. In this disserta¬ 
tion, we also studied the reconstruction of the stau contribution to fc^ by using the information 
which is available at ILC. It was estimated that the contribution can be reconstructed at ~ 0.5% 
at the sample point, which is comparable to or smaller than the measured value of the Higgs 
coupling. Thus, it is possible to test directly whether the excess originates in the stau contri¬ 
bution. Here, the measurement of the stau mixing angle is crucial. We also argued that, if the 
stau mixing angle is measured at the early stage of ILC, it is also possible to predict the heav¬ 
iest stau mass, even when the heaviest stau is not yet discovered at the moment. Therefore, 
the stau contribution to fc^, can be probed not only by discovering the lightest stau, but also by 
studying the stau properties. 

In conclusion, we studied the minimal SUSY models that explain the muon g — 2 discrep¬ 
ancy with the Bino-smuon contributions. It was shown that they are expected to be probed by 
LHC/ILC and LFV/EDM complementarily in future. 
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Appendix A 
Vacuum Decay 


In this appendix, we briefly review the vacuum decay. The vacuum state is stable if the vacuum 
expectation value of scalar potential is at a true minimum. However, when it is at a local 
minimum which is higher than the true minimum, the false vacuum becomes meta-stable and 
could decay into the true vacuum by quantum tunneling effect. The vacuum transition from 
the false vacuum to the true vacuum can be evaluated by semiclassical technique [136,137]. 

The vacuum decay rate comes from the imaginary part of energy of the false vacuum state. 
In general, the energy of the ground state Eq is estimated by path integral method in Euclidean 
space-time as follows. 



(A.1) 


Here, [0] is an Euclidean action, and (f) denotes a scalar field in a theory. In this dissertation, 
(j)i = hy, (j )2 = li, and ^3 = respectively. The integral (A.l) is over all fields satisfying the 
conditions as 



(A.2) 


where, (f)^ is a field configuration of the false vacuum. Eor a stable state, Eq is real. When the 
state is unstable, a non-zero imaginary part appears in Eq, because such unstable state is not 
eigenstate of the Hamiltonian. The decay rate T corresponds to the imaginary part as 


(A.3) 


r = —2 • ImEg- 


The path integral is dominated by bounce configurations, which is a stationary point of the 
action and satisfy the boundary condition (A.2). We assume an 0(4) symmetric action to look 
for the bounce solutions as 


S£[^(r)] = 271^ drr^ 


J 


c 



2 


+n0) 


(A.4) 
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where r is radial coordinate in four-dimensional spacetime, and V((/)) is the potential energy. 
The equation of motion are evaluated as 


3dh„ dV{(f)) 

dr 2 ^ r dr dh^ 

2 ^^+ 

dr^ r dr dcfy^ 


(t = 2,3), 


(A.5) 

(A.6) 


where the boundary conditions are lim^-^o^ = (pi and ^(0) = 0. The solutions in (A.5) 
and (A.6) are the bounce solutions 0;. In this dissertation, they are calculated by numerical 
calculations. 








Appendix B 
Bino Coupling 


In this appendix, we briefly summarize the non-decoupling corrections to the Bino couplings 
with the smuons. The Bino-muon-smuon interactions are 


1 ^ ^ - 
-^int = + h.c., 


(B.l) 


where are the Bino couplings. 

The non-decoupling SUSY-breaking effects are understood in term of renormalization group 
equations. Above the heavy superpartner scale SUSY is exact symmetry. Thus, the rela¬ 
tion gi p^ = gy is maintained. Below the scale the heavy superparticles are decoupled. 

The light particle (SM fermions, light sleptons) loops still renormalize the gauge boson wave 
function. On the other hand, the heavy superparticles loops decouple gauge boson and gaug- 
ino wave function renormalization, respectively. Therefore, SUSY is broken in gauge sector, 
and U(l)y gauge coupling gy and Bino coupling ^ start to evolve differently. 

First, we derive the coefficients Ab in (4.11) and (4.12) of Sec. 4.1.2. At one loop, the 
solution of the RGE of the U(l)y gauge coupling and the Bino coupling are [126] 

1 

gy("tsoft) 

1 


+ 


/3 


gy 


gy(^^soft) 
1 


+ 


8n^ 
Pg 


log 


■log 


-^soft 

5 

^soft 

(B.2) 

-^soft 

5 

^soft 

(B.3) 


where and g are beta function of U(l)y gauge and Bino coupling, respectively. Since 
g^^(MsofJ = gy(Msofj) at boundary, the ratio of (B.2) and (B.3) becomes 


§L,Ri.^soft) 

gyi^soh) 


~ 1 + 


gyC^soft) 

1671^ 




■^soft 

5 

^soft 


(B.4) 


where is the light superparticle scale. Then, let us define 


af — iVsuc2)i ^ 


(B.5) 
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where N^in 2 )i, is a dimension of 811(2)^ representation, N(^ is a color number, and Yj is l/(l)y 
hypercharge of a particle /, respectively. For example, Uq = 2 x 3 x (1/6)^ for the left-handed 
quark doublet. In our setup, the matter content is just SM fermions plus Ugjgpjoji sleptons. The 
beta function of U(l)y gauge and Bino coupling is estimated as^ 

2 r . 

Psy = 3 [3 X (aQ + au + ao + ai + a£) + 2a^J 

If r ^ 1 41 1 

|_^slepton ^e) ^ 

3 

Pg ^slepton y^L ^e) 2 ^slepton" 

Thus the coefficient Ah is evaluated as 

Ab = l3g^- 13^ = ^- (B.8) 

Then, we consider G{a- 2 ) correction due to Wino decoupling. As for logarithmic correction, 
it can be also understood ^(a 2 ) RGEs for the Bino coupling with the left-handed smuon g^. 
If SUSY is exact, the Wino and the W boson contributions to the Bino coupling gi cancel 
each other at 0{a-2). When the Winos are decoupled, only the W boson contribution remains, 
whereas, the U(l)y gauge coupling gy does not run at 0{a-2). Therefore, below the Wino mass 
scale Mgofj, the relation g^ = gy does not hold. Up to G(a 2 ), renormalization group equation 
for the Bino coupling is evaluated as (cf. Ref. [129]) 


(B.6) 

(B.7) 


dg,^ _«2 9 _, I3g 

dlogju 27r4^^ 8n^’ 


(B.9) 


where the first term of right hand side corresponds to the W boson contribution. Solving (B.9) 
and substituting (B.2) for the solution of (B.9), we can obtain Eq.(4.11) in Sec. 4.1.2. 


^ For the full MSSM, coefficients or the beta functions are 

2 p -I 

Pgy = - [3x(aQ-Fau-Fao-Fai-Fa£)-F2aHj 

1 p 

+- [3 X (uq + au + + 2appJ = 11, 

Pg = 3 X {aQ +ay + aj^ + + ag) + 2aH = ll, 

thus, gi n = gy is maintained. 





Appendix C 

Slepton Production 


In this appendix, we briefly summarize the total cross section of slepton production at 
collider. Smuons and Staus are produced via s-channel y/Z exchanges. The production cross 
section is given by [161] 


— 8na^ 3 

cr(e+e = ——Ai 

■' 3s 


■iP_+L^ + 


‘•'sin'^2e^y ^ 

A. 


+ 5;,—(P_+ +P+_) + 5;C; 




-{P.+L + P^_R) 


w 


where we denote £ = fJL, z, and the parameters are defined as 


(C.l) 


A = [1 - (mj^ + mj^f/s] [1 - (mj^ - mj^f/s], (C.2) 

AZ = (C.3) 

s-ml 

Cii /22 = ^[iL+R)±iL-R)cos29j], (C.4) 

1 

Ci 2 = C 2 i = -{L -R)sm26j, (C.5) 

L = -l/2 + sin^0^v, (C.6) 

R = sin^Ow. (C.7) 


The beam polarizations are parameterized as P,,± = (l=FP.-Xl±Prt). 

Selectrons can be produced not only via s-channel y/Z exchanges, but also by the t-channel 
Bino exchange. In this dissertation, we neglect the left-right mixing of the selectrons because 
of small Yukawa coupling. The polarized cross section is evaluated as [161] 


2na^ 


,1 _ ~_ 1 _ ^ ZiVLU. 3 r , -I 2 167rct 

^^^R,L^L,R^ ^L,R^L,R^ ~ ^L,R^z\ I / i / i 
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L,R^R,L ^L,R^L,r) 

2na^ 

3s 

167ra' 


s 

^(^R^R ^L^R^ ~ 

cr(e^e 


['^+gLJigR,L^zY , 


' ^ LA^ra^RB^ lb^ 


V v* v v* 

A=1 B=1 
—> e^~ 1 
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-1 + 2 sin^ 0 , 


gi 

fA 

-^AB 
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sin 0 


2 cos 0 
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w 

2 02 


COS 0iy ’ 


A^-/3^ A^ + ^ 
A^/3-^-^log 


A^-^’ 


-2^i + A^log^ (/i = B) 

fB fA 

H: 


GT = - 


2 


= 


S Aa i Ac 


4^i 


-■AB 


Ag + ^ Aa + 13 

log-T ± log-- , 


(A = B) 
(A/B) 


For the case of Ci rCj^ p^ pairs, the parameter A^ is given by 

2 


4m? 

/3 = 1 -A^. 


= — fm~ — mio^ — 1 , 
5 V ^i.R XV 


while for e^ep^ pairs 

^ = [1 - + m^jVs] [1 - (m^, - A^ = - (^m| + m| - 2 m^o) 

The matrices X^i py^ are given by 


^LA — ^gliON^Al’ — '\/2g^(0^)^i, 


(C.9) 

(C.IO) 

(C.ll) 

(C.12) 

(C.13) 

(C.14) 

(C.15) 

(C.16) 

(C.17) 

1. (C.18) 

(C.19) 


where Ojv is the unitary matrices which dianonalize the neutralino mass matrices. 
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